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Abstract. 

A key source for LISA will be the inspiral of compact objects into massive 
black holes. Recently Mino has shown that in the adiabatic limit, gravitational 
waveforms for these sources can be computed by using for the radiation reaction 
force the gradient of one half the difference between the retarded and advanced 
metric perturbations. Using post-Newtonian expansions, we argue that the 
resulting waveforms should be sufficiently accurate for signal detection with LISA. 
Data-analysis templates will require higher accuracy, going beyond adiabaticity; 
this remains a significant challenge. 

We describe an explicit computational procedure for obtaining waveforms 
based on Mino's result, for the case of a point particle coupled to a scalar field. 
We derive an explicit expression for the time-averaged time derivative of the 
Carter constant, and verify that the expression correctly predicts that circular 
orbits remain circular while evolving under the influence of radiation reaction. 
The derivation uses detailed properties of mode expansions, Green's functions 
and bound geodesic orbits in the Kerr spacetime, which we review in detail. This 
paper is about three quarters review and one quarter new material. The intent is 
to give a complete and self-contained treatment of scalar radiation reaction in the 
Kerr spacetime, in a single unified notation, starting with the Kerr metric, and 
ending with formulae for the time evolution of all three constants of the motion 
that are sufficiently explicit to be used immediately in a numerical code. 
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1. Summary and overview 

1.1. Gravitational radiation reaction of point particles 

A key unsolved problem in general relativity is to compute the gravitational radiation 
produced by a small object spiralling into a much larger black hole. This problem 
is of direct observational relevance. Inspirals of compact objects into intermediate 
mass black holes (M ~ 10 2 — 10 3 Af Q ) may be observed by LIGO and other ground 
based interferometers [I] ; recent observations suggest the existence of black holes in 
this mass range |21 E] • In addition, a key source for the space-based gravitational wave 
detector LISA is the final epoch of inspiral of a stellar-mass compact object into a 
massive (M ~ 10 6 Mq) black hole at the center of a galaxy. Gair et. al. |3] have 
estimated that LISA should see over a thousand such inspiral events during its multi- 
year mission lifetime, based on Monte Carlo simulations of the dynamics of stellar 
cusps by Freitag [H]- 

Observations of these signals will have several major scientific payoffs 6 : 

• From the observed waveform, one can measure the mass and spin of the central 
black hole with fractional accuracies of order 10 -4 [7JIH1- The spin can provide 
useful information about the growth history (mergers versus accretion) of the 
black hole 0. 

• Likewise, one obtains a census of the inspiralling objects' masses with precision 
~ 10~ 4 , teaching us about the stellar mass function and mass segregation in the 
central parsec of galactic nuclei. 

• The measured event rate will teach us about dynamics in the central parsec of 
galaxies. 

• The gravitational waves will be wonderful tools for probing the nature of black 
holes in galactic nuclei, allowing us to observationally map, for the first time, the 
spacetime geometry of a black hole, and providing a high precision test of general 
relativity in the strong field regime ^3 ■ 

To be in the LISA waveband, the mass M of the black hole must be in the range 
10 5 Mq £ M £ 1O 7 M [Tg. The ratio between the mass fi ~ (1 - 100)M Q of the 
compact object and M is therefore in the range 1CP 7 <; fJ,/M <J 10~ 3 . These systems 
spend the last year or so of their lives in the very relativistic regime close to the black 
hole horizon, where the post-Newtonian approximation has completely broken down, 
and emit ~ M/ fj, ~ 10 5 cycles of waveform \12\ . 

Realizing the above science goals will require accurate theoretical models 
(templates) of the gravitational waves. This is because the method of matched filtering 
will be used both to detect the signals buried in the detector noise, and to measure 
the parameters characterizing detected signals. The accuracy requirement is roughly 
that the template should gain or lose no more than ~ 1 cycle of phase compared to 
the true waveform over the ~ 10 5 cycles of inspiral. These sources must therefore 
be modeled with a fractional accuracy ~ 10~ 5 . The past several years have seen 
a significant research effort in the relativity community aimed at providing these 
accurate templates. 

To date, there have been several approaches to generating waveforms. The 
foundation for all the approaches is the fact that, since fi/M <C 1, the field of the 
compact object can be treated as a linear perturbation to the large black hole's 
gravitational field. On short timescales, the compact object moves on a geodesic of 
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the Kerr geometry, characterized by its conserved energy E, z-component of angular 
momentum L z , and Carter constant Q. Over longer timescales, radiation reaction 
causes the parameters E, L z and Q to evolve and the orbit to shrink. 
The various approaches are: 

(i) Use of post-Newtonian methods: Fairly crude waveforms can be obtained using 
post-Newtonian methods |131 |Hj • These have been used to approximately scope 
out LISA's ability to detect inspiral events 0] and to measure the waveform's 
parameters [H]- However, since the orbital speeds are a substantial fraction of 
the speed of light, these waveforms are insufficiently accurate for the eventual 
detection and data analysis of real signals. 

(ii) Use of conservation laws: In this approach |14H15llTClll7| one uses the Teukolsky- 
Sasaki-Nakamura (TSN) formalism ^1^11201 to compute the fluxes of energy E 
and angular momentum L z to infinity and down the black hole horizon generated 
by a a compact object on a geodesic orbit. Imposing global conservation of energy 
and angular momentum, one infers the rates of change of the orbital energy and 
angular momentum. For certain special classes of orbits (circular and equatorial 
orbits) this provides enough information that one can also infer the rate of change 
of the Carter constant Q, and thus the inspiralling trajectory. 

(iii) Direct computation of the self-force: In this more fundamental approach one 
computes the self-force or radiation-reaction force arising from the interaction of 
the compact object with its own gravitational field. A formal expression for this 
force in a general vacuum spacetime in terms of the retarded Green's function was 
computed several years ago Translating this expression into a practical 
computational scheme for Kerr black holes is very difficult and is still in progress. 
Roughly 100 papers devoted to this problem have appeared in the last few years; 
see, for example Poisson [S3] for an overview and references. 

(iv) Time-domain numerical simulations: Another technique is to numerically 
integrate the Teukolsky equation as a 2+1 PDE in the time domain |241 1251 
1261 1271 f2~Sl 1291 130j . and to model the compact object as a finite-sized source. This 
approach faces considerable challenges: (i) There is a separation of timescales 

— the orbital period is much shorter than the radiation reaction timescale. (ii) 
There is a separation of lengthscales — the compact object is much smaller than 
the black hole, (iii) The self-field of the small object must be computed with 
extremely high accuracy, as the piece of the self-field responsible for the self- 
force is a tiny fraction (~ /i/M) of the divergent self-field. This approach may 
eventually be competitive with (ii) and (iii) but is currently somewhat far from 
being competitive [flux accuracies are ~ 10% as compared to ~ 10~ 6 for (ii)]. 

1.2. Radiation reaction in the adiabatic regime 

A key feature of these systems is that they evolve adiabatically: the radiation reaction 
timescale ~ M 2 //i is much longer than the orbital timescale ~ M, by a factor of the 
inverse M/fi of the mass ratio. This has implications for the nature of the signal. The 
self-acceleration a of the compact object can be expanded in powers of the mass ratio 
as 
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Here <Zj diss an d a\ 



,00118 



are the dissipative and conservative pieces of the leading- 



order self-acceleration computed in Refs. Similarly, a 2j diss and a 2iC ons are the 

corresponding pieces of the first correction to the self-acceleration, which has not yet 
been computed (although see Ref. f° r work in this direction). 

The effect of the dissipative pieces of the self force will accumulate secularly, while 
the effect of the conservative pieces will not. Hence the effect of the dissipative pieces 
on the phase of the orbit will be larger than that of the conservative pieces by a factor 
of the number of cycles of inspiral, ~ M / fx. Consider now, for example, the azimuthal 
phase (f)[t) of the orbit. This can be expanded in powers of the mass ratio using a 
two-timescale expansion as l3"3"] 



where the leading-order phase <f>i(t,t\) is determined by ai,di ss , and the higher order 
correction <j>v(t, t\) to the phase is determined by 

0*1 cons cLIld Q-2,diss- 

Here t\ is a "slow" 

time variable which satisfies dt\jdt = 0(fj,/M). 

We shall call leading-order waveforms containing only terms corresponding to 
4>i{t,t\) adiabatic waveforms. From Eq. (|1.2() these waveforms are accurate to ~ 1 
cycle over the inspiral (more precisely, the phase error is independent of fi/M in the 
limit fi/M — > 0). To compute adiabatic waveforms we need only keep the dissipative 
piece of the self-force, and we can discard any conservative pieces. The conservation- 
law method discussed above yields waveforms that are accurate to this leading order, 
but is limited in that it cannot be used for generic orbits. 

Adiabatic waveforms will likely be sufficiently accurate for detecting inspiral 
events with LISA g2|. The correction (f>2 to the phase of the waveform can be 
estimated using the post-Newtonian approximation, and amounts to <J 1 cycle over the 
inspiral for typical parameter values. While the post-Newtonian approximation is not 
strictly valid in the highly relativistic regime near the horizon of interest here, it suffices 
to give some indication of the accuracy of adiabatic waveforms. This computation is 
described in | Appendix A| 

Recently Yasushi Mino derived a key result that paves the way for computations 
of adiabatic waveforms for generic inspirals in Kerr |34) . Mino showed that the time 
average of the self-force formula of [211 for bound orbits in Kerr yields the same 
result as the gradient of one half the difference between the retarded and advanced 
metric perturbations. This "half retarded minus half advanced" prescription is the 
standard result for electromagnetic radiation reaction in flat spacetime due to Dirac 
|35| . This prescription was also posited, without proof, for scalar, electromagnetic 
and gravitational radiation reaction in the Kerr spacetime by Gal'tsov |36| . However, 
prior to Mino's analysis, the prescription was not generally thought to be applicable 
or relevant in Kerr. 

As one might expect, knowledge of the infinite-time-averaged self-force is sufficient 
to compute the leading-order, adiabatic waveforms [HI], although there are subtleties 
related to the choice of gauge [371 EHO EE3 • The sufficiency can also be established 
using a two-timescale expansion |33| . 

1.3. Implementing Mino's prescription for generic inspirals 




) 



(1.2) 



In this paper we apply Mino's result to compute an explicit expression for the time- 
averaged time derivative of the Carter constant. We specialize to the case of a particle 
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endowed with a scalar charge, coupled to a scalar field. This computation is useful as 
a warm up exercise for the more complicated case of a particle emitting gravitational 
waves |%U] . 

We start in Sec.|21by describing the model of a point particle coupled to a scalar 
field, and review how the self-force causes both an acceleration of the particle's motion 
and also an evolution of the renormalized rest mass of the particle |3J. In Sec. 03 we 
review the properties of generic bound geodesic orbits in the Kerr spacetime. A crucial 
result we use later is that the r and 9 motions are periodic when expressed as functions 
of a particular time parameter we call Mino time |34|. and that the t and cj> motions 
consist of linearly growing terms plus terms that are periodic with the period of the 
r- motion, plus terms that are periodic with the period of the ^-motion 02] • 

Section 0] reviews Mino's derivation of the half-retarded-minus-half-advanced 
prescription for radiation reaction in the adiabatic limit [33|. In Sec. we discuss 
a convenient basis of modes for solutions of the scalar wave equation in Kerr, namely 
the "in" , "out" , "up" and "down" modes used by Chrzanowski 03] an d Gal'tsov . 
We also review several key properties of these modes including symmetry relations 
and relations between the various reflection and transmission coefficients that appear 
in the definitions of the modes. Section reviews the standard derivation of the mode 
expansion of the retarded Green's function, and Sec. reviews the derivation of the 
mode expansion of the radiative Green's function given by Gal'tsov correcting 
several typos in Gal'tsov which are detailed in |Appcndix B| 

Next, in Sec. |H1 we turn to the source term in the wave equation for the scalar 
field. We review the derivation of Drasco and Hughes of the harmonic decomposition 
of this source term in terms of a discrete sum over frequencies in which harmonics 
of three different fundamental frequencies occur 02] • We derive expressions for the 
mode coefficients in the expansion of the retarded field near future null infinity and 
near the future event horizon. These coefficients are expressed as integrals over a torus 
in phase space, which is filled ergodically by the geodesic motion. 

Section El combines the results of the preceding sections to derive expressions 
for the time-averaged rates of change of two of the conserved quantities of geodesic 
motion, namely the energy E and the angular momentum L z . The expressions are 
derived from the radiative self force, following Gal'tsov 131]. Gal'tsov also shows 
that identical expressions are obtained by using the fluxes of energy and angular 
momentum to infinity and down the black hole horizon. [This result has recently been 
independently derived for circular, equatorial orbits in Ref. [S|.] We also show that 
the time-averaged rate of change of the renormalized rest mass of the particle vanishes. 

In section ^5] we derive an expression for the time-averaged rate of change of the 
Carter constant, using the radiative self-force and the mode expansion of the radiative 
Green's function. This expression [Eq. Hl().l|) below] is the main new result in this 
paper. It involves two new amplitudes that are computed in terms of integrals over the 
torus in phase space, just as for the amplitudes appearing in the energy and angular 
momentum fluxes. Finally, in Sec. ^2 we show that our result correctly predicts the 
known result that circular orbits remain circular while evolving under the influence of 
radiation reaction. This prediction serves as a check of our result. 

As apparent from the above summary, about 25 percent of this paper is new 
material, and the remaining 75 percent is review. The intent is to give a complete 
and self-contained treatment of scalar radiation reaction in the Kerr spacetime, in a 
single unified notation, starting with the Kerr metric, and ending with formulae for 
the evolution of all three constants of the motion that are sufficiently explicit to be 
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used immediately in a numerical code. 
2. The scalar field model 

We consider a point particle of scalar charge q coupled to a scalar field We denote 
by Ho the bare rest mass of the particle, to be distinguished from a renormalized 
rest mass /i which will occur below. The particle moves in a spacetime with metric 
g a p which is fixed; we neglect the gravitational waves generated by the particle. The 
worldline of the particle is x a = z a (r), where r is proper time. The action is taken to 
be 

S = -\J <* 4 z^(V<f) 2 - Jdr {fi - q<P[z(r)}} . (2.1) 

Varying this action with respect to the worldline yields the equation of motion 

( M o - q$)a a = q(g af3 + u'VjV/j*, (2.2) 

where u a is the 4- velocity and a a is the 4-acceleration. Following Poisson [23], we 
define the renormalized mass fi by| 

H(r) = iia -&[*(?)]. (2.3) 
The equation of motion l|2.2ll can then be written as 

a a = g,a0 + u «j/) V/3 $. (2.4) 
A* ' 

It is also useful to rewrite the definition 1)2. 3[1 of the renormalized mass in terms of a 
differential equation for its evolution with time: 

^(r) = -qu a V a t>. (2.5) 

The phenomenon of evolution of renormalized rest-mass is discussed further in Refs. 
|45l 1461 14*7] . Equations 1)2.4(1 and l|2. 5JI can also be combined to give the expression 

/" = vPVpOiu") = qX7 a $. (2.6) 

for the total self force. The components of f a parallel to and perpendicular to the 
four velocity yield the quantity dfi/dr and /i times the self- acceleration, respectively. 
Varying the action 1)2. 1|) with respect to the field $ gives the equation of motion 

n$(x)=T(x), (2.7) 

where the scalar source is 

/oo 
dTS^[x,z(T)}. (2.8) 

Here 

x — x ')IV—g (2-9) 

is the generalized Dirac delta function, where 5^\x) = S(x°)S(x 1 )5(x 2 )5(x 3 ). We will 
assume that there is no incoming scalar radiation, so that the physical solution of the 
wave equation l|2.7fl is the retarded solution 



*ret(aO = / d A x'^/~g(x')G Ict (x,x')T{x') = -q / drG xet [x, z{t)]. (2.10) 



\ If we endow the particle with an electric charge e and add the electromagnetic coupling term 
e J d,TU a A a to the action 12.11 . the equation of motion 12.41 becomes [ia a = q{g a @ + u a u^)\7 + 
eF a ^up, where F a @ is the Faraday tensor. This shows that the renormalized mass 12.31 is the mass 
that would be measured by coupling to other fields. 
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Here G ret (x,x') is the retarded Green's function for the scalar wave equation 12.711 . 

The retarded field l|2.10|l must of course be regularized before being inserted into 
the expression (|2.4|) for the self- acceleration and into the expression <|2.5|l for d/i/dr, 
or else divergent results will be obtained. The appropriate regularization prescription 
has been derived by Quinn 41 . The regularized self-acceleration at the point z a (r) 
is jH] 

a a (r) = -{g aP + u a uP) \\R M v? - lim f ° dr'V G Tet [z(r), z{t')}} , (2.11) 
and the regularized expression for dfi/dr is 

= -^ R + 1 2 Um f T E dT'u a V a G Ict [z(T),z(T% (2.12) 
dr 12 ^J-oo 

As discussed in the introduction, we specialize in this paper to bound motion 
about a Kerr black hole of mass M with M ^> fi. In this case the terms in Eqs. (I2.11|l 
and (|2.12|l involving the Ricci tensor R a p vanish. Also, as long as the orbit is not 
very close to the innermost stable orbit, the evolution of the orbit is adiabatic: the 
orbital evolution timescale ~ M 2 / '/x is much longer than the orbital timescale ~ M. 
This adiabaticity allows a significant simplification of the formulae l|2.11|l and (|2.12fl 
for the self-acceleration and for the rate of change of mass d/i/dr, as shown by Mino 
|34| . Namely, these formulae reduce to the simple expressions (|2.4() and Q2.5[l . but 
with $ replaced by the radiative field 

$rad = ^($rct-*adv)- (2.13) 

Here $ a dv is the advanced solution of Eq. (|2.7() . In Sec. 0] below we review Mino's 
argument, specialized to the scalar case. 

To summarize, the starting point for our analysis is the self-acceleration expression 

a a = 9 (3 a/3 + u « u /3) V/3 $ rad . (2.14) 

A 4 

together with the expression for the evolution of rest mass 

5r( T ) = -<Z^ Q V Q $ rad . (2.15) 
dr 

Equations l|2.14|) and (|2.15|) are equivalent to the equation 

f a - u p Vp(jiu a ) = gV Q $ rad . (2.16) 
for the total self-force. 



3. Generic bound geodesies in the Kerr spacetime 



This section summarizes the notation and results of Drasco & Hughes 021 > an d a ls° 
some of the results of Mino [31], for generic bound geodesic orbits in Kerr. 
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3.1. The Kerr spacetime 

In Boycr-Lindquist coordinates (t,r,Q,<f)), the Kerr metric is 

, 9 / 2Mr\ l9 4a sin 2 8Mr , ,, , , A , 9 ^sin 2 6> 9 
ds 2 = - 1 — dt 2 dtty + (tn 4 - Aa 2 sin 2 9)— — d(j> 2 

+ Ed6> 2 + -|dr 2 . (3.1) 

Here 

£ = r 2 + a 2 cos 2 0, (3.2) 

A = r 2 + a 2 - 2Mr, (3.3) 

tu = vV +a 2 , (3.4) 

and M, a are the black hole mass and spin parameter. Throughout the rest of this 
paper we use units in which M — 1, for simplicity. The square root of the determinant 
of the metric is 

y/— g = E sin 6, (3-5) 
and the wave operator is given by 



4 

7X7 2-2/1 

— a 2 sm 2 

A 



+ (A$, r ), r + — (sin0$, 9 ), fl . (3.6) 

The differential operator that appears on the right hand side of this equation is the 
Teukolsky differential operator J5j for spin s = 0. 

We will use later the Kinnersley null tetrad I, n, fh, fh* , which is given by 

r=^-d t +d r + ^d ct> , (3.7) 

n = — d t - —d r + —8*, (3.8) 
2£ 2E 2£ 



and 



V / 2(r + zacos 
The corresponding one-forms are 



1 f i 

iasinOdt + de + -r—xd^ I . (3.9) 



and 



1 = -dt + a sin 2 0a> + — dr, (3.10) 

A , aAsin 2 6> 1 , 

n= dt-\ dd> dr, (3.11) 

2£ 2£ V 2 ' 17 



m = — — (— iasin 6dt + Y,de + iw 2 sin Odd) . (3.12) 

V2(r + iacos6) V ' 



The basis vectors obey the orthonormality relations l-ri — —1 and fh ■ fh* = 1 while all 
other inner products vanish. The metric can be written in terms of the basis one-forms 
as 

g a p = -2l^ a n^ + 2m (a m* py (3.13) 
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3. 2. Constants of the motion 

We define the conserved energy per unit rest mass // 

* = --|< (3 ' 14) 
the conserved z-component of angular momentum divided by fj,M 

_ d . 
L z =u-—, (3.15) 

and Carter constant divided by fi 2 M 2 

Q = u 2 e -a 2 cos 2 6E 2 + cot 2 8L 2 Z + a 2 cos 2 6. (3.16) 

[From now on we will for simplicity call these dimensionless quantities "energy" , 
"angular momentum" and "Carter constant" .] The geodesic equations can then be 
written in the form 021 

(J) = [E(r 2 + a 2 ) - aL z ] 2 - A[r 2 + (L z - aE) 2 + Q] = V r (r) , (3.17) 
dO \ 2 

— ] = Q-cot 2 6 L 2 -a 2 cos 2 9{l- E 2 ) = V e {9) , (3.18) 
dX J 

f x = esc 2 6L Z + aE (^±^ - l) - = V<t> (r, 9) , (3.19) 

dt _ 
d\- E 



(r 2 + a 2 ) 2 n ^ 2 t 



— a sin 



aL z [l J—)=V t {r,6). (3.20) 



Here A is the Mino time parameter |34| . related to proper time t by 

dX = -dr. (3.21) 
S 

Also these equations define the potentials V r (r), Vg(9), Vj,(r, 6) and Vt(r,9) §. 

Sometimes it will be convenient to use instead of the Carter constant Q the 
quantity 

K = Q+{L z -aE) 2 . (3.22) 

For convenience we will also call this quantity the "Carter constant". In the 
Schwarzschild limit Q and K are given by Q — L 2 + L 2 and K = L 2 + L 2 + L 2 . 
The quantity K can be written as 

K = K af} u a Uf3, (3.23) 
where K a/3 is the Killing tensor 

K afi = 2Y,m {a m^ - a 2 cos 2 9g af3 . (3.24) 
Using the identity (|3.13l) this can also be written as 

K afi = 2£Z( Q n« + r 2 g afi . (3.25) 

Using the formulae Ij3.7|l and (|3.8|l for the null vectors I and n together with the 
definitions Ij3.14|l and l|3.15|l of E and L z , we obtain from Eq. (|3.25J) the following 
formula for K: 

K = ^(w 2 E-aL z ) 2 - Aw 2 - r 2 . (3.26) 

§ These quantities were denoted R(r), @(0), <E>(r, 8) and T(r,8) in Ref. [421 . We no not use this 
notation here since it would clash with the functions R, and <t> defined in Eq. 15.11 below. 
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Solving this for u r gives 



u 




(3.27) 



this formula will be useful later. 



3. 3. Parameterization of solutions 

Following Mino |34|. we parameterize any geodesic by seven parameters: 
E, L z , Q, A,q, Aeo, to,(j)Q. 



(3.28) 



Here to and 4>o are the values of t(A) and 0(A) at A = 0. The quantity A r o is the value 
of A nearest to A — for which 7'(A) = r m j n , where r m ; n is the minimum value of r 
attained on the geodesic. Similarly Ago is the value of A nearest to A = for which 
9(A) = 9 min , where 9 m i n is the minimum value of 9 attained on the geodesic. This 
parameterization is degenerate because of the freedom to reparametrize the geodesic 
via A — > A + AA. We discuss this degeneracy further in Sec. 13.71 

Frequently in this paper we will focus on the fiducial geodesic associated with the 
constants E, L z and Q, namely the geodesic with 



3.4- Motions in r and 9 

It follows from the geodesic equations l|3.17|l and (|3.18|) that the functions r(A) and 
9(A) are periodic. We denote the periods by A r and Ag, respectively, so 



Using the initial condition r(A r o) = ?~min we can write the solution r(A) to Eq. (|3.17|l 
explicitly as 



Aj-o — -^eo — to — <f>a — 0. 



(3.29) 



r(A + A r ) =r(A), 



9(X + Ag) = 9(X). 



(3.30) 



r(A) =f(A-A r0 ), 
where the function f(A) is defined by 




(3.31) 



(3.32) 



Similarly we can write the solution 6(A) to Eq. I|3.18|l explicitly as 

0(A) = §(A - Aeo), 
where the function 9(A) is defined by 



(3.33) 




(3.34) 



The functions f(A) and 9(A) are just the r and 9 motions for the fiducial geodesic. 
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3.5. Motion in t 

Next, the function V t (r, 9) that appears on the right hand side of Eq. (|3.2U|) is a sum 
of a function of r and a function of 9: 

Vt(r,e) = V tr (r) + V te (6), (3.35) 

where V tr (r) = Em 4 / A + aL z (l - zu 2 /A) and V te (9) = -a 2 Esm 2 9. Therefore using 
t(0) = to we obtain 

t(X) = t + [ d\' {V tr [r(X')] + V te [9(X')}} . (3.36) 



Next we define the averaged value {Vt r ) of Vj r to be 



i r A " i ^ 



(V tr ) = -p / dX V tr [r(X)\ = — dX V tr [f(X)]. (3.37) 

"■r Jo "r Jo 

Here the second equality follows from the representation H3.31() of the r motion together 
with the periodicity condition l|3.30|l . Similarly we define 



(Vte) = -T- / dX V te [9(X)] = — / dX V te [9(X)]. (3.38) 
A<? Jo A e Jo 

Inserting these definitions into Eq. (|3.36(l allows us to write t(X) as a sum of a linear 

term and terms that are periodic 42 : 

t(X) =t a + TX + At r (X) + At e (X) (3.39) 

= t + TX + At(X). (3.40) 
Here we have defined the constant 

T = (Vtr) + (Vte), (3.41) 
and the functions 

At r (X) = ^ dX'(vtrlr(X')}-(Vt r )\, (3.42) 

Atg(X) = J dX'{v te [e(X')} - (V te )\. (3.43) 

The key property of these functions is that they are periodic: 

At r {X + A r ) = At r (X), (3.44) 

At e (X + Ag) = At fl (A); (3.45) 

this follows from the definitions 13.42JI and H3.43|l together with the periodicity 
condition (|3.3U|) . We can exhibit the dependence of these functions on the parameters 
X r0 and Ago by substituting the formulae l|3.31|l and (|3.33|) for r(X) and 9(X) into Eqs. 
lEHU and The result is 

Atr(A) = i r (X - X r0 ) - ir(-X r o), (3.46) 

Atg(X) = t e (X ~ X ea ) - t e (-X eo ), (3.47) 

where the functions t r (X) and tg(X) are defined by 

t r (X)= [ dX'(v t r{r(X')]-(Vt r )\, (3.48) 



r-A 

tg(X)= I dX'{Vte[9(X')]- (V te )\. (3.49) 
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3.6. Motion in 

The motion in can be analyzed in exactly the same way as the motion in t. First, 
the function V$ (r, 9) that appears on the right hand side of Eq. (|3.19|) is a sum of a 
function of r and a function of 9: 

V4r,9) = V^ r (r)+V^(0), (3.50) 
where V,p r (r) = aE(w 2 /A — 1) — a 2 L z /A and V$g(6) = csc 2 8L Z . Therefore using 
0(0) = 4>o w e obtain 

0(A) = O + f dX' {^ r [r(A')] + V<pg[9(\')]} . (3.51) 
Jo 

Next we define the averaged value (V0,,) of to be 

(V*) = T" / dA ^V[r(A)] = 1T f ' dX W(A)]. (3.52) 

Jo "-r Jo 

Here the second equality follows from the representation (|3.31|) of the r motion together 
with the periodicity condition l|3.30|l . Similarly we define 



(V^) = -T- / dA V* fl [0(A)] = — / dA V>[0(A)]. (3.53) 
Ae Jo A e Jo 

Inserting these into Eq. I|3.51[) and using 0(0) = 0o allows us to write 0(A) as a sum 
of a linear term and terms that are periodic |42j : 

0(A) =0 o + T A + A0 r (A) + A0 e (A) (3.54) 

= 0o + T A + A0(A). (3.55) 

Here we have defined the constant 

t* = 0V) + (M> (3.56) 

and the functions 

A0 r (A) = ^ dA'{vV[r(AO]-(VV>}, (3.57) 

A0#(A) = £ dX'(v^[9(X')} - (V»j. (3.58) 

The key property of these functions is that they are periodic: 

A0 r (A + A r ) = A0 r (A), (3.59) 

A0 e (A + A e ) = A0 e (A); (3.60) 

this follows from the definitions 13.57fl and l|3.58|l together with the periodicity 

condition (|3~30"|) . We can exhibit the dependence of these functions on the parameters 

A r o and A#o by substituting the formulae l|3.31(l and l|3.33(l for r(A) and 9(X) into Eqs. 
(|T5?jl and The result is 

A0 r (A) = r (A - A r0 ) - ^r(-Aro), (3.61) 

A0e(A) = 9 (A - A eo ) - <M~A eo ), (3-62) 

where the functions r and 0# are defined by 

$ r (X)= [ dX'\v$ r [r(X')] - (V$ r ) j, (3.63) 



e (A)=/ dX'\v^[9(X')]-(V^)\. (3.64) 
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3. 7. Re-parameterization freedom 

Not all of the parameters E, L Z1 Q, \ rQ , Aeo, t and <pQ that characterize the geodesic 
are independent. This is because of the freedom to change the dependent variable 
A via A — * A = A + AA. Under this change of variable the parameters A r o and Aeo 
transform as 

A r o — * A r o = A r o + AA, (3.65) 

Aeo — * A eo = Aeo + AA. (3.66) 

We can compute how the parameters to and 4>o transform as follows. Combining Eqs. 
H3.40fl . (|3.46|l and l|3.47|l gives the following formula for the t motion: 

t = t + TX + t,,(X - X r0 ) - i r (-X r0 ) + t e (X - Aeo) - te(-*oo)- (3-67) 

Rewriting the right hand side in terms of A, X r0 and Aeo yields 

t = t + TX- TAX + t r (X - X r0 ) - i r (-X r0 ) + *e(A - X eo ) - ig(-X eo ) 

= t + TX + t r (X - X r0 ) - i r (-X r o) + i e {X - Ago) - 4(-A eo ). (3-68) 

Comparing the first and second lines here allows us to read off the value of to: 

t = h- TAX + i r (-X r0 - AX) - 4(-A r0 ) + i {-Xg O - AA) - i e (-X 80 ). (3.69) 

Similarly we obtain 

0o -> 4>o = 0o-T AA + I .(-A rO -AA)-0 r (-A r o) + 0e(-Aeo-AA)-00(-Aeo).(3.7O) 

In Sec. 18. 41 below we will explicitly show that all of the amplitudes and waveforms we 
compute are invariant under these transformations. 

3.8. Averages over Mino time in terms of angular variables 

We will often encounter functions Fg(X) of Mino time A that are periodic with period 
Ag. For such functions the average over Mino time is given by 

(Fg)x = ~t~ [ " F e (X)dX. (3.71) 
A e Jo 

In order to compute such averages it will be convenient to change the variable of 
integration from A to a new variable x- This variable is a generalization of the 
Newtonian true anomaly and is defined as follows The potential Vg(8) defined 

by Eq. (j3.18|l can be written in terms of the variable 

z = cos 2 <9 (3.72) 

as 

1 — z 

Here j3 = a 2 (l — E 2 ) and z_ and z + are the two zeros of Vg(z). They are ordered such 
that < z_ < 1 < z + . The variable \ is defined by 

cos 2 §(X) = z_ cos 2 x, (3.74) 

together with the requirement that x increases monotonically as A increases. From 
the definition l|3.34|l of 6(X) together with the formula (|3.73|) for Vg(z) we get 

^ = ^P(z+-z- cos 2 x). (3.75) 
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Therefore the average (|3.71|) can be written as 

{F e )x = -r- / d X ■ (3.76) 

A similar analysis applies to averages of functions F r (X) that are periodic with 
period A r gj. The average over A of such a function is 



(F r )x = "t / F r (X)d\. (3.77) 

A r Jo 

Now the potential V r (r) in Eq. 13.17|l can be written as 021 

V r {r) = (1 - E 2 )( n - r)(r - r 2 ){r - r 3 )(r - r 4 ), (3.78) 

where r%, r 2 , r% and r± are the four roots of the quartic V r (r) = 0, ordered such that 
T/k < t% < r% < rj.. For stable orbits the motion takes place in r2 < r < r\. The 
orbital eccentricity e and semi-latus rectum p are defined by 

n = ^ (3.79) 



and 



r-2 = j?-. (3.80) 

1 + £ 



We also define the parameters p% and p^ by 

r* = — e , ^ = — (3.81) 
Then following Ref. ^1] we define the parameter ip by 

f(A) = V - -. (3.82) 

W 1 + £ COS ij) ^ ' 

It follows from these definitions that 02] 

^=W), (3-83) 

where 



V(ip) = ^} f 2 b(l-£)- P3 (l + £cosV')] 1/2 [p(l + £)~P4(l+£cos^] 1/2 . (3.84) 

1 — £^ " 

The average (I3.77JI over A is therefore 

(, r>l . > _ Jff-WW). (3. 85) 

Note that the viewpoint on the new variables x,"0 adopted here is slightly different 
to that in 021 : Here they are defined in terms of the fiducial motions f(A) and 
Q{\) instead of the actual motions r(A) and 0(X). The new viewpoint facilitates the 
computation of amplitudes and fluxes for arbitrary geodesies; if one is working with 
the fiducial geodesic the distinction is unimportant. 
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3.8.1. Bi-periodic functions Finally suppose we have a function F(X r ,Xg) of two 
parameters A r and Xg which is biperiodic: 

F(X r + A r , Xg) = F(X r , Xg) 

F(X r ,Xg + Ag) — F(X r , Xg). (3.86) 
The average value of this function is 



< F >A = ^r/ dX r / 0lXgF(X r ,Xg). (3.87) 

By combining the results l|3.76[) and l|3.85[) we can write this average as a double 
integral over the new variables x an( i V*" 



21 F[X r (^),Xg(x)} 



VP( Z + z ~ cos2 X)F{ip) 
This formula will be used in later sections. 



(3.88) 



4. Mino's derivation of half retarded minus half advanced prescription 

In this section we review the proof by Mino |34j that for computing radiation reaction 
in the adiabatic limit in Kerr, one can use the "half retarded minus half advanced" 
prescription. We specialize to the scalar case. 

We start by defining some notation for self- forces. Suppose that we have a particle 
with scalar charge q at a point V in a spacetime M with metric g a p. Suppose that 
the 4- velocity of the particle at V is u a , and that we are given a solution <I> (not 
necessarily the retarded solution) of the wave equation l|2.7|l for the scalar field, for 
which the source is a delta function on the geodesic determined by V and u a . The 
self- force on the particle is then some functional of V, u a , g a f} and $, which we write 
as 

f a [P,u a ,g al3 ,$] . (4.1) 

Here we suppress the trivial dependence on q. Note that this functional does not 
depend on a choice of time orientation for the manifold, and also it is invariant under 
u a -u a . 

Next, we define the retarded self-force as 

fZt [V, u a , g a p] = f a [P, u a ,g a p, $ rct ] , (4.2) 

where $ re t is the retarded solution to the wave equation (|2.7|l using the time orientation 
that is determined by demanding that u a be future directed. Similarly we define the 
advanced self force by 

ia"dv [P,u a ,g aP ] = r [P,u a ,g aP ,<t> adv } , (4.3) 

where $ a dv is the advanced solution. It follows from these definitions that 

fZt [V, -u a ,g a0 ] = /« dv [P, u a ,g a p] . (4.4) 
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4-1- Properties of the self- force 

The derivation of the half-retarded-minus-half-advanced prescription in the adiabatic 
limit rests on two properties of the self-force. The first property is the fact that 
the self-force can be computed by subtracting from the divergent field <E> a locally 
constructed singular field $ s i ng that depends only on the spacetime geometry in the 
vicinity of the point V and on the four velocity at V . This property is implicit in the 
work of Quinn and Quinn and Wald [22], and was proved explicitly by Detweiler 
and Whiting 05]. We can write this property as 

f a [V, u a , g a p, $] = qV a ($ - $ Bing [P, u a , g af3 }) . (4.5) 

The second property is a covariance property. Suppose that tp is a diffeomorphism 
from the manifold M to another manifold iV that takes V to ^("P). We denote by i/j* 
the natural mapping of tensors over the tangent space at V to tensors over the tangent 
space at tp(V). We also denote by ip* the associated natural mapping of tensor fields 
on M to tensor fields on N. Then the covariance property is 

fW7') ) f/,f W ,f$] = ff[?' 1 «», W) $], (4.6) 

This expresses the fact that the self-force does not depend on quantities such as a 
choice of coordinates. It follows from the definition (|4.2[) that the retarded self-force 
/ r " t satisfies a similar covariance relation: 

ftWJ'l.Wfw] = rf" e t[P,u a ,g af} }. (4.7) 



4-2. Property of generic bound geodesies in Kerr 

Next we review a property of bound geodesies in Kerr upon which the proof depends. 
This is the fact that for generic bound geodesies there exists isometries ip of the Kerr 
spacetime of the form 

t^2t 1 -t, (p -> 2fa - fa (4.8) 

where t\ and fa are constants, which come arbitrarily close to mapping the geodesic 
onto itself. 

To see this, note that if a geodesic is mapped onto itself by the mapping (|4.8() . 
then the point [t, r(t), 9(t), <f)(t)] is mapped onto [2ti — t, r(t), 8(t), 2fa — fat)], which 
must equal [t' , r(t') , 9(t') , fat')] where t' = 2t% — t. If we specialize to t = t\ then we 
see that fat\) = fa. We denote this point by Q — (ti,n,9i,fa). Next, since the 
geodesic is determined by Q and by the four-velocity u a at Q, the geodesic is mapped 
onto itself by ip if and only if ip*u a — —u a at Q. This will be true if and only if the 
t and 9 components of u a vanish. Thus, a geodesic will be invariant under a map of 
the form (j4.8(l if and only if it contains a point Q which is simultaneously a turning 
point of the r and 9 motions. 

Consider now a generic bound geodesic. Such a geodesic can be characterized by 
the parameters E, L z , Q, A r o, A#o, to> 4>o m the notation of Sec. 13.31 The turning 
points of the r motion occur at A = X r .„ = A r0 + nA r /2, where n is an integer. 
Similarly the turning points of the 9 motion occur at A = Ag TO = Xg + mAg/2, where 
m is an integer. Since generically A r and Ag will be incommensurate, for any e > 
we can find values of the integers m and n so that 

|A nn - A<?, m | < e. (4.9) 
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Thus we can find a point on the geodesic that is arbitrarily to close to being a turning 
point for both the r and 9 motions, and hence there is an isometry of the form l|4.8|l 
that comes arbitrarily close to mapping the geodesic onto itself. 

Below we will simply assume that the isometry does map the geodesic onto itself; 
the error associated with this assumption can be made arbitrarily small in the adiabatic 
limit 1 1 . 



4-3. Property of the adiabatic limit 

The self-acceleration a a at a point ?ona geodesic has three independent components, 
since a a u a — 0. These three components are determined by the time derivatives of 
the three conserved quantities 

dE dL z dK 
dt ' dt ' dt 

In order to compute the evolution of the orbit in the adiabatic limit, it is sufficient to 
know the time averaged quantities 

'dE\ /dL z \ / dK \ 

~dt/t \~te/t \te/t 
where 

dE\ 1 f T , dE . . , A . 

— ) = m — / dt—(t). 4.11) 

dt It T-.00 2T J_ T dt y ' y ' 

This was shown in Ref. ,34 (although see Refs. [^TJESIESI f° r caveats related to gauge 
issues), and is also discussed in detail using a two-timescale expansion in Ref. 



(4.10) 



4-4- Derivation 

We now combine the various results discussed above to derive the half retarded minus 
half advanced prescription. We parameterize the geodesic as x a — z a (t) , and denote 

by 

f I % t (t) = fr a et[z a (t),U a (t) >9a0 } (4.12) 

the retarded self force at z a (t). The key idea is to use the covariance property l|4.7|l . 
with t/j taken to be the isometry associated with the geodesic. If we take V to be 
z a (t), then the left hand side of (|4.7II evaluates to 

f« t [z a (2h - 0,-" a (2ti - t),g aP ] - / a Q dv [^(2ti - t),u a (2h - t),g a0 ] 

= f^{2h-t). (4.13) 

Here we have used the property l|4.4|) and the fact that ip*g a (3 = g a p. This gives 

/^ v (2t 1 -t)=^/r a et(*)- (4-14) 

Now let £ denote one of the three conserved quantities E, L z or K. The time 
derivative of £ depends linearly on the self force, so we can write 

§(*)=H a (t)/£ t (i) (4.15) 

j Note however that real inspirals will consist of finite curves which are approximately locally 
geodesic; the corresponding error for such curves will contribute to the post-adiabatic correction 
to the inspiral, i.e., to the second term in Eq. 11.21 . 
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for some B a . Below we will compute S Q explicitly for the three different conserved 
quantities and derive the key property that 

1>*E a (t) = -E a (2t 1 -t). (4.16) 
Now acting on Eq. 1(4.15(1 with ip* and using the relations (|4.14(l and 14.16(1 yields 

(Z a fZt)(t) = -(H a /« dv )(2t 1 - t). (4.17) 
Taking a time average gives 

( S a/ret)t = ~ (^aiTdv)* ' (4-18) 

and hence 

§) = (S Q / r Q ct) t = ~ (2 Q (/ r Q ct - /a Q dv)> t • (4-19) 

Finally we use the explicit formula ((4.5(1 for the self force. The contributions from the 
locally constructed singular field <t> s i ng cancel, and we obtain 

'^^|(s Q v Q (a> rct -a> adv )) t . (4.20) 



This can be written as 
'd£\ 
dt/ t 

where 



(Ha/£d) t (4.21) 



/rad = |v a ($ I . t-$ a dv) (4.22) 

is the self-force one obtains from the half retarded minus half advanced prescription, 
cf. Eqs. lEHHl and (|2~TK|l above. What we have shown is that although the quantity 
/ r " d differs from the true self- force / r " t , the difference averages to zero when one takes 
a time average over the entire geodesic. 

It remains to derive the formula ((4.161) . For energy and angular momentum, we 
can write the conserved quantity as the inner product of a Killing vector with the 
4-velocity of the particle: 

£ = i a u a . (4.23) 

Here £ = —d/dt for £ — E, and £ = d/d<j> for £ — L z . Taking a derivative with 
respect to proper time r gives 

— = u^V^u^) = {Vpi a )u f3 u a + Uu Vpu a ) = H a a a , (4.24) 

where a a is the 4-acceleration. Here we have used the fact that £ is a Killing vector so 
V( ct £ / 3) = 0. Next we use the geodesic equations 1(3. 2UH and 1(3.21(1 to translate from a 
proper time derivative to coordinate time derivative, and we use the relation between 
self- force and self-acceleration given by Eqs. ((2.4(1 and 1(2.6(1 . Using the definition ((4.15(1 
of S a now gives 

The formula l(4~Tol) now follows from ip*u a (t) = -u Q (2ii - t) and ip*£ a = -£ a . 
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Turn now to the Carter constant. Taking a time derivative of the expression 
CEIIII for K gives 

^ = u'V^iKapu"^) = V h K a/3) u^u a u f) + 2K a pu a v?V~ f u p 

= 2K af} u a a 13 . (4.26) 

Here we have used the Killing tensor equation V( 7 if Q( g) = 0. Proceeding as before 
now gives the expression 

2Ti(r) 

Ea = Ja^T) {S « + u ^) K ^ ul - ( 4 - 27 ) 

The formula (|4.16|) now follows from the relations tp*u a (t) = ~u a (2ti — t) and 

ip*K al3 = K at3 . 

5. Separation of variables and basis of modes 

5.1. Separation of variables for scalar wave equation 

Separation of variables for the scalar wave equation in Kerr was first carried out by 
Dieter Brill and others in 1972 |5l)|. and subsequently generalized to higher spins by 
Teukolsky ^J^l- We substitute the ansatz 

r, 6, <j>) = R(r)e{6)e m ' t 'e^ luJt (5.1) 

into the homogeneous version of Eq. 12. 7f) . and make use of the expression (|3.6(l for 
the wave operator. The result is the two equations 

— !— (sin ee e )e- -^5-6 + a 2 uj 2 cos 2 6*9 = -A9, (5.2) 
smo ' " sin 

and 

(AR r ) tr + — (lo 2 w a - Aramuj + m 2 a 2 ) R = (A + ui 2 a 2 )R. (5.3) 

Here A is the separation constant We label the successive eigenvalues A of Eq. 1)5. 2f> 
by an integer I, with I = \m\, \m\ + 1, \m\ + 2, . . ., and we denote these eigenvalues by 
A w / m - In the Schwarzschild limit a = we have A w ; m = 1(1 + 1). Since the differential 
equation (|5.2|l does not depend on the sign of to or the sign of m, we have 

A- w ;- m = A w ; m , (5-4) 

a property which will be used later. We denote by 9 w / m (#) the corresponding solutions 
to Eq. (|5.2|) . The functions 

s ulm {0A) = e im<t >e Urn (e) (5.5) 

are the spheroidal harmonics. These harmonics are orthogonal on the sphere, and we 
can choose them to be orthonormal like the spherical harmonics: 



d Q S u im{9, 4>)*S u i> m >(6, 4>) =5wS mm i. (5.6) 
Following Gal'tsov we also choose the phases of the spheroidal harmonics to satisfy 

(PS ulm )(9,(f>) EE S u l m (w- 0,7T + <£) = (-l) l S ulm (e,4>) (5.7) 

and 

S- ut - m (0, 4>T = {-l) m S w i m (e, (f>). (5.8) 
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Here P is the parity operator. 

The radial equation (|5.3|l can now be simplified by defining 

u{r) = wR{r), (5.9) 
and by using the tortoise coordinate r* defined by 

dr* = (vu 2 /A)dr. (5.10) 
An explicit expression for r* is 

r*= r + 2r+ In T -=I± - ^ In (5.11) 



where r± — 1 ± \J\ 2 — a? are the two roots of A(r) = 0. [Recall that we are using 
units in which the mass M of the black hole is unity.] The resulting simplified radial 
equation is 

d 2 u 

^{r*) + V ulm [r*)u{r*) = Q, (5.12) 
where the potential is 

A ,„ o as 3r 2 A 2 2 
V u im = =■ 3r 2 - 4r + a 2 + =- + u? 

H j [to 2 a 2 - Aramuj ~ A(A w ; m + uj 2 a 2 )] . (5.13) 

From Eq. (|5.4|l it follows that this potential satisfies 
u)l— m 

(r) = K J / m (r), (5-14) 
an identity which will be used later. 

5. 2. Basis of modes 

In this section we follow the treatment of Chrzanowski 03] as modified slightly by 
Gal'tsov . We follow Gal'tsov's notation except that we omit his subscripts s 
denoting spin, since we have specialized to s = 0. 

We start by analyzing the asymptotic behavior of the potential (|5.13|) as r* — > — oo 
and as r* — > oo. We consider first the limit r* — > — oo, the past and future 
event horizons. In this limit A — > and r — > ?*+ = 1 + \/l — a 2 . We find 
VLim(r) Vui m (r+) = pL. where 

Prow = U3 — TOW+, (5.15) 

and 

W+ = JL (5.16) 
2r + 

is the angular velocity of rotation of the horizon. Therefore the solutions of the radial 
equation near the horizons are of the form 

u{r) oc e ±iPm " r \ (5.17) 

In the limit of r* — > oo (past and future null infinity), VLi m (r*) — > cj 2 , so the solutions 
are of the form 

u(r) oce ±ia,r *. (5.18) 



Computing inspirals in Kerr in the adiabatic regime 



21 



5.2.1. "In", "up", "out" and "down" modes We now define, following Gal'tsov, the 
following solution 

( T i It) |-l/ 2 p-iPmur* r * _> _ no 

This equation serves to define the mode as well as the complex transmission and 
reflection coefficients r w / m and o^im ■ The coefficient a u i m is a normalization constant 
whose value is arbitrary; we will discuss a. convenient choice of OLujim 

later. We will 

often denote the set of subscripts uilm collectively by A: thus o~\, t\, and sometimes 
we will omit the subscript and use simply a, r. 

The "in" mode (|5.19|l is a mixture of outgoing and ingoing components at past 
and future null infinity, since the mode function is multiplied by e~ lu)t . At the past 
and future event horizons, the mode is purely ingoing when the sign of p mul is the 
same as the sign of uj. However, from the definition (|5.15fl of p mul we see that the 
sign of ujpmuj can be negative; this occurs for superradiant modes. Thus, at the future 
event horizon the "in" modes can be either ingoing or outgoing. 

The crucial feature of the "in" modes is that they vanish on the past event horizon. 
This feature will be used later in constructing the various Green's functions. A more 
precise statement of the result is that a solution $ of the wave equation which is 
a linear combination of "in" modes with coefficients c w ; m , such that the coefficients 
depend smoothly on u) (a reasonable requirement), must vanish at the past event 
horizon. To see this, note from Eqs. (|5.1|1 and l|5.9|) that the solution can be written 
as 

[■ca 00 ' in ( *\ 

$(t,r,e,ct>)= / dioJ2 Y, e ' luJtc ^S^ m (0,^) ""7 V ■ (5.20) 

We now insert the asymptotic form l|5.19|l of the mode function near the horizon, and 
we use the definition l|5.15|l of p muJ . This gives 

1 f°° 

${t,r,6,(t>) = - / duj Ve-^( t+I '*)c (i)/ „ l S'^ m (0,0)a w;m T tj;m | J5mw r 1 / 2 e lm ^ r *(5.21) 
V^ + r*;^™^ 1 -*. (5.22) 

uj 

Ira 

Now all of the quantities that depend on lu in the integrand are smooth functions of 
ui. Since Fourier transforms of smooth functions go to zero at infinity, it follows that 
the function Gi m (v]Q,4>) defined by Eq. I|5.22|l satisfies G; m — > as v — > — oo, where 
v = t + r* . Thus, will vanish as v — > — oo, on the past event horizon. 
Another important property of the "in" modes is the identity 

<W-m(r*)=^(rT- (5-23) 

To derive this, note that the potential V w / m is real and also satisfies the identity 
H5.14JI . It follows that the function u™ ul _ m (r*)* satisfies the same differential equation 
H5.12JI as u™ lm (r), and has the same asymptotic behavior as r* — > oo if wc choose the 
normalization constants to satisfy 

ot-u-m = o£i m . (5.24) 
The result l|5.23|l now follows. In particular the identity (|5.23(l implies that 

(5.25) 



tm 

i x ^ 
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Figure 1. An illustration of the various types of modes in black hole spacetimes. 
Here J~ denotes past null infinity, future null infinity, E~ the past event 
horizon, and -B+ the future event horizon. The four panels give the behavior of the 
four different modes "in", "out", "up" or "down" as indicated. A zero indicates 
the mode vanishes at the indicated boundary. Two arrows indicates that the 
mode consists of a mixture of ingoing and outgoing radiation at that boundary. 
Two arrows with an "OR" means that the mode is either purely ingoing or purely 
outgoing at that boundary, depending on the relative sign of p m[J and lo. The 
"in" modes vanish on the past event horizon, and the "up" modes vanish on past 
null infinity. Thus the "in" and "up" modes together form a complete basis of 
modes. Similarly the "down" and "out" modes together form a complete basis of 
modes. 



and 



T—wl—m — T ujlm- 

We next define the "up" modes 



(5.26) 



up 



lm( r ) — P^ lr - 



-1/2 up,, 



7 P 
I III ( 



-ip-n 



- DO 



(5.27) 



As before this equation defines both the mode as well as the complex coefficients \x^i m 
and Vojim- The coefficient (3 u i m is a normalization constant whose value is arbitrary; 
we will discuss a convenient choice of /3 w / m later. The "up" modes are a mixture of 
ingoing and outgoing components at the past and future event horizons. At future 
null infinity, the mode is purely outgoing. By using an argument similar to that given 
above, one can show that these modes vanish at past null infinity. Also we can argue 
as before that 



j I in 



( r *y, 



(5.28) 
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M-a>!-m = /Cim! ( 5 - 29 ) 

and 

V-ul-m = Kim' ( 5 - 30 ) 
as long as we choose the normalization constants to satisfy 

P-ull-m = fiiohn- (5.31) 

Since the "in" modes vanish at the past event horizon, and the "up" modes vanish 
at past null infinity, the two sets of modes are orthogonal to one another and together 
form a complete basis of modes. 

Next we note that (u™ lm )* is a solution of the differential equation (|5.12|) since 
the potential H5.13fl is real. We denote this quantity by 

<tm = (<l m y- (5-32) 
Similarly we define 

utZ n = K P iJ*- (5-33) 
These "down" and "out" modes together form a basis, as the "out" modes vanish on 
the future event horizon and the "down" modes vanish at future null infinity. The 
properties of all of these modes are summarized in Fig. ^ 

5.2.2. Relations between scattering and transmission coefficients The four sets of 
constants o^im, r w / m , \i^im and v u lm are n °t all independent. One can derive relations 
between these coefficients by using the fact that the Wronskian 

W{U 1 ,U 2 ) = "l-TT - ~^J U2 - ( 5 - 34 ) 

is conserved for any two solutions u\ and u 2 of the homogeneous equation Ij5.12|l . By 
evaluating W(u^ m , u"\ m ) at r* = — oo and at r* = oo using the asymptotic relations 
(|5.19[) and l|5.27|l and equating the results, we obtain 

Tujlm^ujlm = 1- (5.35) 

A similar calculation with the modes w u ,?„, and yields 

T ulm v ^lm — ~°~ulm- (5.36) 

Finally, using the modes u™ lm and yields 

|^ m | 2 + ^^|r. ;m | 2 = l. (5.37) 

5.2.3. Relations between mode functions The "down" and "out" modes can be 
expressed as linear combinations of the the "in" and "up" modes, since the latter 
form a basis of modes. Using the asymptotic forms Ij5.19|l and (|5.27|l of the modes at 
r* — > oo together with the definition (|5.33|l allows us to identity the coefficients for 
utZ n , giving 

u dow„ (r) = Sdin u ^ m{r) _ aulm ^hn u ^ m(r) . (5 . 38) 

A similar computation using the definition (|5.32(l and using the asymptotic forms at 
r* — > — oo gives 

•7/ out f<r\ — ^P™^ T ^lm a uilm „.up t \ T ulm V ulm "ulm ^in ( \ (r o Q \ 
u u)lm\ r ) — 1 1 -3 U uilm\ r ) u uilm\ r )-\°- iiy ) 

This can be simplified using the relations (|5.35(l and H5.36|l to give 



W = r^^l^P^LW + ^<im< lm (r). (5.40) 
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5.2.4- Complete mode functions Finally, following Gal'tsov |36| . we define the 
complete mode functions 

2 



7r "L(*> r > 6, >< 



2 

— t 

ZU 

2 

— ( 

w 







(5.41) 






(5.42) 


Sulm{9, 




(5.43) 



and 



n d JZ n (t, r, 9, 0) = -e^S^O, 0)t$T(r*). (5.44) 
w 

Gal'tsov actually includes a factor of 1 +pP in these definitions, where P is the parity 
operator defined by 

(Pf)(t, r, 9, 0) = f(t, r,ir-9,n + ( f>) (5.45) 

for any function / and p = ±1 is an additional mode index. However when acting on 
the scalar modes considered in this paper, we have from Eq. I|5.7J) that 

l+pP = l+p(-l) l = 25 p<i _ 1) i. (5.46) 

Thus the p index is redundant in the scalar case; we have simply dropped it and 
replaced the factor of 1 + pP by 2. 

6. Construction of the retarded Green's function 

6.1. Formulae for retarded Green's function 

The retarded Green's function G ret (x,x') is defined such that if $ obeys the scalar 
wave equation with source T 

□$=T, (6.1) 

then the retarded solution is 



$ret(z) = / d 4 xW-9(x')G Tet (x,x')T(x'). (6.2) 



The expression for the retarded Green's function in terms of the modes defined in the 
last section is% 



-°° 1=0 m=-l a " l ™P" 1 ™ 

x [<ln^Xf m ^rO{r - r>) + C Jm (s)7$T (zTflCr' - r)] . (6.3) 

Here 9(x) is the step function, defined to be +1 for x > and otherwise. Note 
that the expression (|6.3|) is independent of the values chosen for the normalization 
constants a u i m and /3^/ m , since the factor of 1/a cancels a factor of a present in the 
definition (|5.19|) of the "in" modes, and similarly for j3 and the "up" modes. 

% This expression agrees with Eq. (2.17) of Gal'tsov l-jlil . if we assume that the normalization 
constants obey oiuilmPujlm = ~ 1/(16tt). However Eq. (2.25) of Gal'tsov says ctujlmPujlm = 1/4; thus 
we disagree with Gal'tsov by a factor of — Air. For this reason we do not assume any values for the 
normalization constants a^im and /3 w ( m in our computations, but leave them unevaluated. 
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The expression (|6.3[1 can be expanded into a more explicit form by using the 
definitions (|5.41|) - (|5.44|) of the complete mode functions 7r w ; m (t, r, 9, 0) in terms of 
the radial mode functions M w ; m (r), together with the definitions i|5.32[l and l|5.33|l of 
the "out" and "down" modes. This gives 

G rct (t,r,8,<P;t',r',9',<b') = — / fa-^'V^fMl^^^Tn 

x ; 1 n K p lm (r)<i m (r')6(r r>) + < lm {r)uJ lm {r')9{r' - r)] . (6.4) 

In this section we review the standard derivation of the formula l|6.4JI . 

6.2. Derivation 

The key idea behind the derivation is the following. Suppose that the source T(x) is 
non-zero only in the finite range of values of r 

(6.5) 

Then, the retarded solution $ rct (a;) will be a solution of the homogeneous equation in 
the regions r < r mul and r > r max . Now, the retarded solution is determined uniquely 
by the condition that it vanish on the past event horizon and on past null infinity. 
This property will be guaranteed if: 

(i) When we expand $ re t in the region r < r min on the basis of solutions ^ lm (x) 
and 7r"^„ of the homogeneous equation, only the "in" modes contribute. Then, 
since the "in" modes vanish on the past event horizon, $ rot must also vanish on 
the past event horizon. 

(ii) When we expand <& rc t in the region r > r max on the basis of solutions Tr™ lm (x) 
and 7T^? m , only the "up" modes contribute. Then, since the "up" modes vanish 
on past null infinity, $ rc t must also vanish on past null infinity. 

We start by defining the Fourier transformed quantities 



T(«, r, 9, (/)) = I dte wt T(t,r,e,(t)) (6.6) 

and 

r, 9, <t>) = I dte wt <S>(t,r,6,(t)). (6.7) 



[For convenience we drop the subscript "ret" on <f> in the remainder of this section.] 
We make the following ansatz for the Green's function: 

/OO 7 
_ e -ic( t -i')G ret (r-, 9, <f>; r' , 9', </>'; w). (6.8) 

Inserting these definitions into the defining relation (|6.2|l and using the formula Q3.5jl 
for y/— g gives 

r, 0,4)= f dr' J d 2 ft'£(r', 9') G rct (r, 9, 0; /, 9', 0'; w)f (w, /, 9', <f>'). (6.9) 

Next, we decompose the quantities and ET on the basis of spheroidal 
harmonics: 

$(w,r,0,<fi = ^S u i m (6,<t>)Ru,i m (r) (6.10) 
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and 

im(r) (6-11) 

im 

We include the factor of E on the left hand side of Eq. I|6.11[l because of the appearance 
of a factor of E in the integrand of Eq. i|6.9[l . The factor of r 2 on the right hand side of 
Eq. (|6.11l) is so that the coefficients T ulm reduce to the conventional spherical harmonic 
coefficients in the Schwarzschild limit a = 0. From the orthogonality relation (|5.6[l the 
inverse transformations are 



Rvlm(r) = / d 2 n S ulm (9, <jy)* r, 6, <jy) (6.12) 

and 

r 2 %i m {r) = ( d 2 fl S u im(0, <t>T E(r, 6)f(u, r, 8 : 0). (6.13) 



Next, we insert the decompositions (|6.1U|I and (|6.11l) of $ and T into the Fourier 
transform of the differential equation l|2.7|l using Eq. I|3.6I) . This gives 



d u^ijn 



+ Vui m {r*)u u i m (r*) = s w i m , (6-14) 



dr* 2 

where u u i m {r) — n7 (r)i? w ; m (r), the potential V^zm is given by Eq. (|5.13() . and the 
source term is 

r 2 A - 

Sculm = ^Tujlm- (6.15) 

We denote by G ti; / m (r*, r*') the relevant Green's function for the differential equation 

p OO 

dr*'G ulm (r*,r*')s wlm (r*'). (6.16) 

OO 

We will derive an explicit formula for G w ; m (r* , r* ') shortly. First, we note that 
the Fourier-transformed retarded Green's function G re t(r, 0, (f>; r' , 8', <//; oj) can be 
expressed in terms of G^im via the formula 

G rct (r, 0, <f>; r', 6', 0'; u) = V S ulm (9, <P)S^ lm (9', 0')* G ^ r *f\ (6.17) 

* — ' Turn 

lm 

To see this, insert the ansatz (I6.17f) into the relation l|6.9|l and simplify using the 
definition (|6~T^|) of f alm . This gives 

roc ' 2 

Comparing this with the definition (|6.1UI) of R u i m and simplifying using the relations 
<f5~§|) . lf5~TU|l and (|6~T3|) finally yields the formula l|6.16|) . Hence the ansatz (|6~T71 is 
correct. 

Finally we turn to the derivation of the formula for the Green's function G^im 
for the differential equation (|6.14(l . From the discussion at the start of this section, 
the relevant boundary conditions to impose are that 

G ulm {r*,r*')^^ lm {r*), r* -» -co (6.19) 

and 

G u i m {r*,r*')KuZ m {r*), r* -» oo. (6.20) 
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Consider now the expression 
Guim{r*,r*') ' 



+ < lm {r)uZJr')e(r-r') 



(6.21) 



where W is the conserved Wronskian l|5.34|l . This expression satisfies the boundary 
conditions (|6.19|) and (|6.2U|) . Also one can directly verify that it satisfies the differential 
equation 1)6. 14[) with the source replaced by S(r* — r*'), using the fact that the "in" 
and "up" modes satisfy the homogeneous version of the differential equation. This 
establishes the formula 1|6.21[1 . 

Next, we compute the Wronskian W{u™ lm , u^ m ) using the asymptotic 
expressions (I5.19|) and H5.27J) for the mode functions for r* — > oo. This gives 

W(tft lm ,u% m ) = 2ia^ lm ^ lm f- (6.22) 

\0J\ 

Inserting this into Eq. I|6.21|) and then into Eqs. I|6.17|l and l|6.8|l finally yields the 
formula (|6.4|) . 

7. Construction of the radiative Green's function 

7.1. Formulae for radiative Green's function 

Using the retarded Green's function G rc t(x, x') discussed in the last section we can 
construct the retarded solution $ ret (x) of the wave equation l|2.7|) . Similarly, we 
can construct the advanced solution $ a dv(^) using the advanced Green's function 
G a dv(^) x'); this is the unique solution which vanishes on the future event horizon and 
at future null infinity. One half the retarded solution minus one half the advanced 
solution gives the radiative solution: 

$rad (x) = i [$ ret (x) - $ ad v (x)] . (7.1) 

Clearly the radiative solution is given in terms of a radiative Green's function 

$ rad (z) = / d 4 xW-9{x')G ta , d (x,x')T{x'), (7.2) 

where 

G rad (x, x') = i [G rct (x, x') - G adv (x, x')] . (7.3) 

The expression for the radiative Green's function in terms of the modes defined in 
Sec. is 



1 f°° °° LU 
J -°° l=o m=-l 1 1 



UPmu \ T uilm\ 2 down/ \ down r/\* 
>~ulm \ x )~u}lm \ x ) 



OLuln 

(7.4) 



m I 

Note that this expression actually independent of the values chosen for the 
normalization constants a w ; m and (3 u i m , since the factor of l/|a| 2 cancels factors 
of a present in the definition l|5.19|l and l|5.32|l of the "out" modes, and similarly for 
(3 and the "down" modes. 
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The expression <|7.4[1 can be expanded into a more explicit form by using the 
definitions (|5.41|) - (|5.44|) of the complete mode functions 7r w ; m (t, r, 9, 4>) in terms of 
the radial mode functions u^i^r). This gives 

G iad (t,r,8,fct',r',9'^') = — / due^*-*)YS ulm (8,il>)S ulm (0',<j>'"^ 



. i / ^uji in \ ^ i ry ^ujirii it) i / 

87TI /_„„ WUJ 1 
Ira 



1 „out M-,/ 011 * [VI* I l T ojim| down/ \ down/ /\* 



(7.5) 



In this section we review the derivation of the formula (|7.5[1 given in Gal'tsov jS,- 
7.K. Derivation 

We start by deriving the identity 

G aAv {x,x') = G xe t{x',x). (7.6) 
For any two functions $(x) and ^(x) we have the identity 

V Q ($V a * - *V a $) = (7.7) 
Integrating this identity over a spacetime region 1/ and using Stokes theorem gives 

/ ($V a * - *V a $)d 3 S a = / d 4 a;V^g($n* - *□$)• (7.8) 

Jav Jv 
If we choose V to be the region of spacetime outside the black hole, then the boundary 
dV of V consists of the past and future event horizons, and also past and future null 
infinity. 

We now fix two points x\ and x 2 in V, and we choose 

$(x) = G nt (x,xi) (7.9) 

and 

=G, dv (x,x 2 ). (7.10) 
Then □$ = S^ix, x x ) and = 5^(x, x 2 ) , where S^' is the covariant delta function 
H2.9(l . The right hand side of Eq. Ij7.8jl therefore evaluates to 

$(x 2 ) - *Oi) = G Tet (x 2 ,xi) - G adv (xx,x 2 ). (7.11) 
Consider now the left hand side of Eq. I|7.8[) . Since <i> vanishes on the past event 
horizon and at past null infinity, and VP vanishes on the future event horizon and 
at future null infinity, there is no portion of dV where both $ and iff are nonzero. 
Therefore the left hand side vanishes, so the expression H7.11I) must vanish, proving 
the identity (Tfy 

We now insert the identity (|7.6(l into the definition l|7.3(l of G ra d, and use the 
explicit expansion (j6.4|l of G rc t- We also specialize to 

Ctujlm — Pulm = 1; (7-12) 

this causes no loss in generality since the final result (|7.5|> is independent of the values 
of and Puim- This gives 

G iad (x,x') =— diue-^-^Ys^ue^^s^ie'^r— — 7 

8m J_ 00 \w\ wvj' 

x K P i m (r)<im(r')e(r r') + u% hn (r)u% m (r'W - r)\ 



8m J_„ f-^ 



x K P i m (r')<l m (rW -r)+ <UO<LM% - r')} , (7.13) 
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where x = (t,r,9,<f)) and x' = (t' ,r' ,9' ,(/)'). In the second sum we relabel to — * — u 
and m — > — m, and we use the formula (|5.8|) . This gives 

1 /-oo 1 

8ttj „/_ 00 cj row' 



x ( K P i m (r>Zm(r') + uZi-Jr^^ir')] 9(r r') 

+ [< lm {r)nZ m {r') + ^_ m (r)u^_ m (r')] 9{r> -r)\. (7.14) 



Next, since our choice (J7.12JI of the normalization constants satisfies the conditions 
l|5.24[) and l|5.31|) . we can use the identities (|5.23|l and l|5.28|l . This yields 

1 f°o 1 



\U!\ VOW 

Im 



+ Kim(rK P lm (r') + <i m (rTuZ m (ry] 9(r' r) . (7.15) 



Consider now the coefficient of 9(r — r') inside the curly brackets in Eq. (|7.15|) . 
We denote this quantity by H (r, r'): 

H{r, r') = u up (r)u in (r') + u up (r)*u in {r')* . (7.16) 

Here and in the next few paragraphs we omit for simplicity the subscripts ulm. Next 
we use the formula 15.40f) for the "out" modes in terms of the "in" and "up" modes, 
together with our choice l|7.12|l of normalization constants, to obtain 

u up (r) = y^r K ut (r) - aV n (r)l . (7.17) 

| T |2 

Here 

k = ojp mu /\ojp mu \ (7-18) 

is a variable which is +1 for normal modes and —1 for superradiant modes. Inserting 
this into Eq. I|7.16|) and expanding yields 



u ln (r)*u m (r') + u m (r)u m {r')* - o*u ln {r)u ln (r') 



au ln (r)*u ln (r')* 



(7.19) 



Next, by using Eq. I|7.17[) to evaluate the quantity u np (r)*u up (r') we obtain the identity 

|r| 4 w up (r)*w up (r') - u in {r)u in {r')* - \a\ 2 u in (r)* u in (r' ') = -a*u in {r)u in {r') 

-o-u in (r)*u in {r')*. (7.20) 

We now use this identity to eliminate the last two terms inside the square brackets in 
Eq. (|7.19() . and simplify using the identity (|5.37() . This gives 

H(r, r') = u in {r)*u in (r') + K\T\ 2 u up {r)*u np {r'). (7.21) 

Now the right hand side of Eq. (|7.21|) is explicitly invariant under the combined 
transformations of interchanging r and r' and taking the complex conjugate. However, 
from the definition (|7.16(l of H(r,r'), the left hand side is invariant under complex 
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conjugation. It follows that both sides of Eq. I|7.21[l are real and also symmetric under 
interchange of r and r' : 

H(r,r')=H(r',r). (7.22) 

Next, the quantity inside the curly brackets in the expression l|7.15|l for G ra( j is 

H(r, r')9{r - r') + H(r', r)6{r' - r). (7.23) 

Using the symmetry property (|7.22|l together with 0{r — r') + 6(r' — r) = 1, this can be 
written simply as H{r, r'). Therefore we can replace the expression in curly brackets 
in (|7.15() with the expression (|7.21() for H(r,r'). This gives 

1 f°° uj 1 

G rad (x,x') =^-J ^e-^*- t )^S w/m (M)^ m (0',^r ^ 



lm 



u^i m {r)*u^ lm {r') + r— — \T^ m \ 2 ul P lm (r)*ul P lm (r') 

I ^fmw | 



(7.24) 



Finally, rewriting this in terms of the "down" and "out" modes using the definitions 
H5.32J) and (|5.33(l . and reinserting the appropriate factors of the normalization 
constants yields the formula (|7.5(l . 

8. Harmonic decomposition of the scalar source 

8.1. The retarded and radiative fields 

Using the expression l|tj.3[) for the retarded Green's function together with the integral 
expression (|6.2|l . we can compute the retarded field & lct (x) generated by the source 
T{x). For the case we are interested in, T(x) will be nonzero only in a finite range of 
values of r of the form 

r m m < r < r max . (8.1) 

For r > r max only the first term in the square brackets in Eq. I|tj.3f) will contribute, 
and the function 9{r — r') will always be 1. This gives 

i poo oo l 

U-)^//EE Si J-«» r>r m , (8.2) 

where the coefficients Z°^ n are given by the inner product of the "out" modes with 
the source: 

Z Zlm = (nZllli'T) ■ ( 8 - 3 ) 

Here following Gal'tsov 1351 we have defined the inner product of any functions Q(x) 
and ty(x) on spacetime by 



($,*) = J d 4 xy/-g(x)$(x)*y(x). (8.4) 
Similarly for r < r min we obtain 

, „oo oo i 

= — / dc]T E m R— Z ^»irn(x), r < r min , (8.5) 

where 

Z^ = (^Z n ,T)- (8-6) 
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Finally, the expression l|7.4|) for the radiative Green's function together with the 
definitions i|8.3[l and i|8.6[l give the following expression for the radiative field 



2 poo 00 ' 

$ d (a;) — / duj'S^ 

latH ; 32niJ_ 00 ^, 



7°ut out / \ 



^Pmuj \Tujlm\ 2 ^down down f \ 
\UPmu\ \Pulm\ 

All of these expressions depend on the amplitudes and Z^J°™ 



•7) 



8. 2. Harmonic decomposition of amplitudes 



For the case considered here where the source T(x) is a point particle on a bound 
geodesic orbit, the amplitudes and Z^°™ n can be expressed as discrete sums over 
delta functions |4*2*1IKT] : 

00 00 

Z out, down \ \ ^ ^out,down c/ \ /o o\ 

u,lm = 2s 2s Z lmkn 5{W-Uj mkn ). (8.8) 



k= — 00 n— — 00 



Here 



u m kn — rnQ,^ + kftg + nfl r , (8-9) 

fi* = ^*, Sio = jr, 0- = ^, (8.10) 

re = ^, T. r = ^, (8.11) 
Ag A r 

and r and are defined by Eqs. I|3.41|) and (|3.56[1 . The formula for the coefficients 

^out „ J ^down : 
^Imkn dllu ^Imkn lh 

y out4own _ 47Tg -im<t> a iu mkn tg f ' ji / "//L p^TeAe+nT^r) 

x e [e( ^ )re _^A^(A r ) e -^.(A.) 

x e ^ mfen A tr (A r ) e i Wmfe „Ai e (A e ) u out,djwaj r(Ar) j*_ (g 12) 

Here the functions At r , Ai e , A^ r and A</> 9 are defined in Eqs. J2321>, (G33), <|3~5?|) . 
and (|3.58|) . Also the function QuiimiO) is defined by the decomposition of the spheroidal 
harmonic S^im,: 

S u i m (e,ct>) = Qu,i m (d)e im< l', (8.13) 

cf. Eq. I|5.5|l above. For evaluating the double integral (|8.12|) it is convenient to use 
the identity l|3~H5)l that expresses the integral in terms of the variables x an( i V'- 

In this section we review the derivation of the harmonic decomposition l|8.8[) 
given by Drasco and Hughes |42II51| . adapted to the scalar case, and generalized from 
fiducial geodesies to arbitrary bound geodesies. We consider only the "out" amplitude 
Z°^; the "down" case is exactly analogous. 
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8.3. Derivation 

We start by inserting the definition l|2.8[) of the source T(x) into the definition of 
given by Eqs. (|8.3|) and 18.4fl . This gives an expression consisting of an integral along 
the geodesic of the mode function: 

z°JL = -q r dT*2L[t(j),r(T),e(j),<KT)]*. (s.w) 

J — oo 

Inserting the expression (|5.43[) for the mode function K™^(t,r,9,(j)) in terms of the 
radial mode function (r), and using the decomposition 1)8. 13[) of the spheroidal 
harmonic S^imi®, 4>) now gives 

f°° 1 

We next change the variable of integration from proper time r to Mino time A using 
Eq. i|3.21[) . and we use the expressions (|3.39[1 and 13.54fl for the functions t(X) and 
4>(X). This gives 

/>oo 

ZZTm = ~2qe- ml ^e^ t0 / dAe lA ( r "- mT *)S[r(A),0(A)]e w/m [0(A)]* 



x e -imA^(A) e -imA^(A) e ^At r (A) e ^At fl (A) 1 out Lmi^ (gjg) 

nr[r(A)J 

We now define the function of two variables 
•WAr, A e) = - 2ge- im *>e*"° £[r(A r ), 0(A e )]e w/m [0(A e )]*e-"" A ^> 

x e -»mA<» 9 (A 9 ) e ^At r (A r ) e ^At 9 (A 9 ) 1 [r(A r )]*. (8.17) 

This function has two key properties. First, when evaluated at A r = Xg = A, it 
coincides with the integrand of Eq. I|8.16|l . up to the first exponential factor: 

/•oo 

Z°Jrn= dAJ^ m (A,A)e iA ( rw -" lT ^. (8.18) 



Second, the function J w ; m is biperiodic: 

Julm{X r + A r , Xg) = Jujlm{X r , Xg), 

Julm{X r ,Xg + Ag) = J w / m (A r , Xg). (8.19) 

This follows from the fact that the functions r(A), At r (X) and A(p r (X) are all periodic 
with period A r , and the functions 0(A), Atg(X) and A0g(A) are all periodic with period 
Ag, cf. Sec.|5]above. Standard properties of biperiodic functions now imply that J w ; m 
can be written in terms of a double Fourier series: 



oo oo 



■WA r ,A 9 )= J2 E J^mkne^ {kTeXe+nT ^\ (8.20) 

fc= — oo n= — oo 

where Tg and T r are given by Eq. (|8.11|) . and where the coefficients J u i m kn are given 

by 

Ju,*mfe« = -rV / r dX r f S dXge^ kr ^ +nr ^J ulm {X r ,Xe).(8.21) 

A r Ag J J Q 
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We now insert the Fourier series l|8.2U[l for J w / m , evaluated at X r = Xg = A, into 
the expression (|8.18|l for Z°^. This gives 



/oo 
dAe iA(r W -mT^-kT 9 -nT r )j ( 
. . . -oo 



ujhnkn 



2tt 



— ^ — S(w — LUrakn) Jujlmkn, (8.22) 

where we have used the definition 18.9fl of uj m kn- Now comparing with Eq. I|8.8[) we 
can read off that 

out 27T 

Z°mkn = —Jio mkn lmkn- (8.23) 

Combining this with the definitions (|8.17|l and (|8.21|l evaluated at u> — w m k n finally 
gives the formula (|8.12|) . 

Note that it follows from the harmonic decomposition (|8.8|l that for geodesic 
sources, the continuous frequency lo and the discrete indices I, m are replaced with 
the four discrete indices k,n,l, and m. In this context the operation 

lo — ► — u), m — > —to, I — ► Z (8.24) 

associated with the symmetries ljo"3|) . lUTgj l. (UHlf) . lf!T2"3|) - lfS~2l)j) and lf5~2"g|) - (|d1HT) 

is replaced by the operation 

k — > —k, n — ► — n, to — > — m, ^ — > i. (8.25) 
<§.^. Dependence of amplitudes on parameters of geodesic 

The amplitude Z°^ kn is a function of the parameters characterizing the geodesic, 
namely E, L z , Q, to, tfio, X r o, and Ago, cf. Sec. 13.31 above. We write this dependence 
as 

Z°mkn = Z?mkn(E> L Z ,Q, to, <t>0, A r0 , Ago). (8.26) 

We now specialize to the fiducial geodesic associated with the constants E, L z and Q, 
i.e. the geodesic for which to — 4>o = A r o = Ago = 0. For this case we can simplify 
the formula (|8.12() by setting to and 4>o to zero, by replacing the motions r(A) and 
0(A) with the fiducial motions f(A) and (9(A) defined by Eqs. I|3.32|l and l|3.34|l . and 
by replacing the functions At r , Atg, A0 r and A(f>g with the functions t r , tg, <p r , and 
4>g. defined by Eqs. (|3^9"|) . EHSJ, and filSty . This yields 

Z^ kn {E,L z ,Q,Q,Q,Q,Q) = -^L-J\x r j\x e e i(*T^+»r r A r ) 

S[f(A r ),0(Ag)] Q ~ X 1 *.-im#r(Ar) 

w[r(X r )\ 

For more general geodesies, the amplitude Z°^ kn depends on the parameters to, c/>o, 
X r o and Xgo only through an overall phase: 

ZZ\ n (E, L z , Q, t , O , A r0 , Xg Q ) = e *x lmk Ato,4>o,Ko,Xeo) z out k j Ej L ^ Q> Qj Q) Q) Q)) (8 28) 
where 

Xlmkn(to,4>0, A r0 , Xg ) = kTgXg + nT r X r0 + TO r (— A r o) + (t>g( — \eo) ~ 4>0 

— U m kn [tr(-Ko) + tg(~Xg ) ~ t Q ] ■ (8.29) 
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This formula can be derived by substituting the expressions (|3.4(j[l , l|3.47|l , l|3.61|l and 
H3.62fl for the functions At r , Atg, A(f) r and A<pQ into Eq. I|8.12[l . making the changes 
of variables in the integral 

A r — > X r = A r — A r o, Xg — ► Ae = \e — Xeo, (8.30) 

and comparing with Eq. I|8.27[l . Finally we note that the phase (|8.29|) and amplitude 
(|8.28|l are invariant under the transformations i|3.65[l . (|3.66[) . (|3.69[1 . and (|3.7U|) that 
correspond to the re-parameterization A — > A + AA. This invariance serves as a 
consistency check of the formulae, since we expect the invariance on physical grounds. 



9. Expressions for the time derivatives of energy, angular momentum, 
and rest mass 



9.1. Time averages 

Let £ be one of the three conserved quantities of geodesic motion, E, L z or Q. For 
the purpose of evolving the orbit we would like to compute the quantity 

'*>.■ 

that is, the average with respect to the Boyer-Lindquist time coordinate t of the 
derivative of £ with respect to t. However, the quantity that is most naturally 
computed is the derivative with respect to proper time r, and the type of average 
that is most easily computed is the average with respect to Mino time A. In this 
section we therefore rewrite the quantity (|9.1|) in terms of a Mino-time average of 
d£/dr. 

In the adiabatic limit, we can choose a time interval At which is long compared 
to the orbital timescales but short compared to the radiation reaction timc + . Then, 
to a good approximation we have 

'§>,-£ 

where A£ is the change in £ over this interval. Now let AA be the change in Mino 
time over the interval. From Eq. H3.40[l we have 

At = TAX + oscillatory terms. (9.3) 

Now the oscillatory terms will be bounded as At is taken larger and larger, and 
therefore in the adiabatic limit they will give a negligible fractional correction to At. 
Hence we get 

'd£\ _ 1 A£ 
~dt/ t ~ f AA 



~f\dX/ } 



(9.4) 

A 

where the A subscript on the angular brackets means an average with respect to A. 
Note that using the definition l|3.41|l of T we can rewrite this formula as 

d£\ (d£/d\) x 

dt/ t (dt/dX) x ' 1 ' ' 

+ A natural choice for At is the geometric mean of the orbital time and the radiation reaction time; 
this is the time it takes for the phase difference between the geodesic orbit and the true orbit to 
become of order unity. 
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Finally we can use Eq. I|3.21|l to rewrite the Mino-time derivative in Eq. I|9.4|l in terms 
of a proper time derivative. This gives the final formula which we will use: 

'*\ = l -Uf-) . (9.6) 



dt 1 1 r \ dr I 



A 



9. 2. Formulae for time derivatives of energy and angular momentum 

We now specialize to the cases of energy and angular momentum, £ = E or £ = L z 
For these cases the result for the time-averaged time derivative is 



/ — \ — _ / L ° mkn X u)mkn 



\dt/ t 64tt 2 ^^{ m J \ui mkn \ 



1 i/70ut|2 i ^mknPmkn \ t a\ i^down|2 



l^out \Z _|_ tii^iir- iil^/i \^7Q 



\ua\ 2 A \^mknPmkn\ \Pa\ 2 A 

(9.7) 



Here 

Pmkn = Prnu mkn = Umkn ~ TOW+, (9.8) 

cf. Eq. H5.15fl above. Also the symbol A is a shorthand for the set of indices Imkn; we 
define 

«A = Oilmkn = <Xu> mhn lm> (9-9) 

and similarly for []\ and ta- The sum over A is defined as 

oo oo oo I 

E = E E EE • 

A k= — oo n= — co 1=0 m=— I 

The expression in curly brackets in Eq. I|9.7|l means either a factor of u m kn (for energy) 
or a factor of m (for angular momentum). The first term in the square brackets 
corresponds to the flux through future null infinity, and the second term corresponds 
the flux through the horizon. Note that for energy the future null infinity term is 
always negative; this makes sense since the expression (|9.7|) is the rate of change of 
orbital energy. 

The flux expression l|9.7|l has been derived from the radiation reaction force. As 
shown by Gal'tsov, the same result is obtained by computing the flux directly using 
the stress-energy tensor of the scalar field . 

Note that the final answer l|9.7|l is independent of our choices for the normalization 
constants a a and /3a, since the amplitudes Z^ and Z^ own contain factors of a a and 
Pa that compensate for the factors that appear explicitly in Eq. (|9.7|l . 



9.3. Derivation 

For energy and angular momentum, we can write the conserved quantity as the inner 
product of a Killing vector £ Q with the 4-velocity of the particle: 

£ = £ a ii a . (9.11) 

Here £ = —d/dt for £ = E, and £ = d/d(f> for £ = L z . Taking a time derivative gives 

^ = rfV^taU") = (Vp&Ju'V + U^Vpv?) = Si a a a , (9.12) 

where a a is the 4-acceleration. Here we have used the fact that £ is a Killing vector 

so V (Q £/3) = 0. 
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Next, we use the formula l|2.14|) for the self-acceleration. This gives 

f*£ = £ [cv a ^ aA + (&u t Vv /3 $ rBd ] 

rfr \i L 

= H$- + ^. (9.13) 

Now the second term here is a total time derivative to leading order in q, since £ is 
conserved to zeroth order in q. Hence the change in £ over an interval from t± to t-i 
associated with the second term will be 

-£[*«dfo)-*rad(n)]. (9-14) 

A* 

This is a term which will oscillate but will not grow secularly with time (since $ ra d 
does not grow secularly with time). Hence in the adiabatic limit this term will be 
smaller than the change due to the first term in Eq. (|9.13ll by the ratio of the orbital 
timescale to the inspiral timescale, which is negligible in the approximation we are 
using. Dropping the second term in Eq. (|9.13|) and substituting into Eq. i|9.6fl we get 

'd£\ q_ 
dt 1 1 fiT 

Next, we can obtain an explicit expression for the radiative field $ r ad by 
substituting the harmonic decomposition (|8.8II into Eq. (|8.7f) . This gives 

oo OO OO / 

= db £ £ EES 



:srv Q $ r ad) A . (9.i5) 



oa| 2 



fc= — oo 7l=— oo i=0 m=—l 

WmknPmkn | t a| 7 down down / \ 
\2^A n A \ x ) 



(9.16) 



I |j8a| 

Here as above A denotes the set of indices Imkn, and we have defined 

out, down out, down /q * ^\ 

We now substitute this expression into Eq. (|9.15|) . and use the fact that the operator 
£"V Q gives a factor of ito (for the energy) or im (for the angular momentum). We get 



\~dt/ t ~ 32^f Y I m J |Wmfcn| 



Zout /v 1 out\ 
■> A X^^A /; 

OA 



WmknPmkn | t a| ^down /v _down\ 



(9.18) 



||/3a| 2 a x a /a . 
where the expression in curly brackets means either a factor of u> m kn or a factor of m. 
We now discuss the evaluation of the quantity 

<E*ZSLi m > A (9-19) 
which appears in Eq. (|9.18|) . This is the time average of the "out" mode 
function evaluated on the geodesic, and can be easily evaluated using the harmonic 
decomposition of Sec. |S| By comparing the integrands of Eqs. I|8.14() and (|8.18() we 
find that 

nz a (X)KZ n im[^W] = ~-Ju mknlm {\\ye-^— ~ r * m \ (9.20) 
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where z a {\) = [t(X),r(X),9(X),(p(X)] is the geodesic. Now using the Fourier series 
expansion 18. 21) [I of J u i m gives 

£[*°(A)]7rS^ nI Jz"(A)] =~E ■C™»n«m*'n'e- <Ar(w "*»- w '-*'-' ) . (9.21) 

9 fc',n' 

From this expression it is clear that averaging over A kills all the terms in the sum 
except k! — k, n' — n, and we obtain 

(S[z«(A)]^ ;m [z Q (A)]) A = -lj* mhnlmkn = -l-(Z^ kn )*. (9.22) 

For the last equation we have used Eq. (|8.23() . Finally substituting this result, together 
with a similar equation for the "down" modes, in to Eq. (|9.18l) we obtain the final 
result C£3 



9.4- Time derivative of renormalized rest mass 

Equation (|2.15|) for the time derivative of the renormalized rest mass /i(t) can be 
immediately integrated with respect to r between two times n and r 2 . This yields 

m(t 2 ) - n{n) = q^d[z(T 2 )} - g^rad^(ri)], (9.23) 

where x a — z q (t) is the geodesic. Now since the zeroth-order orbit is bound, the right 
hand side of this equation does not contain any secularly growing terms. Instead it 
consists only of oscillatory terms. Therefore, over long timescales and in the adiabatic 
limit (the regime in which we are working in this paper), the rest mass is conserved: 

^ = 0. (9.24) 

GIT 

10. Expression for the time derivative of the Carter constant 

10.1. Formula for time derivative 

We now turn to the corresponding computation for the Carter constant K defined by 
Eq. ipH5)l . The result is 



dK \ _ 1 LO mk 

A 



dt I \ 327T 2 /! \u m kn\ 



1 ( ^out\* ryout , ^mknPmkn | 7"A 1 2 / ydown \ * ry&, 
• \ Z A ) + I" 1 \Ta 12 \ A ) X A 



own 



OAp \LOmknPmkn\ \f3\\ 2 

(10.1) 

This expression has the same structure as the expression (|9.7f) for the time derivatives 
of the energy and angular momentum, except that the squared amplitudes |^X ut | 2 an d 
|^down|2 haye been rep i acec i witri pro ducts of the amplitudes Z^ ut and Z^ own with two 
new amplitudes and Z^ own . These new amplitudes are defined by the equation 

ZZln™ = -^-e-^e*"-*-* / ' dX r [ ' dX e e^ kT ^ +nT ^e^ mknlm [e(^)}* 



4irq 

— t 

e 

— imA0 r (A r ) — imAtf>g(\e) iw mt „At r (A r ) iw m(o „ Atg(Xg) 



TA r A 



o Jo 

f dR ou ^ ,( ^ own 

x G mfc „[A r , A e ]^; d X n [r(A r )]* + G[A r , A e ] [r(X r )}* 



(10.2) 
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where 



out, down/ \ 
u <*JLim ( r ) 



n out,down/ \ w m fc n im v / ,-,,■> , 

Kj knlm W = g(r) , (10-3) 



G mfc „(A,,, Afl) = - — (m 2 E - aL z ){w 2 uj mkn - am) + 2iru r A, (10.4) 

G(A r , Ag) = — iAEu r . (10.5) 

On the right hand sides of Eqs. (|10.4[) and (|10.5|) . it is understood that £(r, 0), A(r) 
and ro(r) are evaluated at r = r(A r ) and 9 = 9(Xg). Also it is understood that u r is 
given as a function of r by Eq. 1)3. 27|) . and hence as a function of A r by using r = r(X r ) 
and by resolving the sign ambiguity in Eq. Ij3.27|l using Eq. I|3.82|l . 

Several features of the formulae I|1U.1|I and (|10.2|) are worth noting: 

• The formula for the ne w am plitudes Z°^ kn and Zf%£. has a very 
similar structure to the formula l|8.12|l defining the amplitudes Z°^ kn and Zf£^. 
In particular, the derivation of the formula (|8.28|) for the dependence of the 
amplitudes on the parameters to, </>o, A r o and Xeo of the geodesic (via an overall 
phase) carries through as before. This means that the time derivative of the 
Carter constant is independent of these parameters, as expected, since the phase 
from the untilded amplitudes cancels the phase from the tilded amplitudes in Eq. 

itnm . 

• Although we have no general definition for "flux of Carter" , presumably the first 
term in Eq. (|10.1fl corresponds to something like the flux of Carter to future null 
infinity, and the second term to the flux of Carter down the black hole horizon. 

• As was the case for the original amplitudes, to evaluate the double integral (|10.2fl 
it is convenient to use the identity (|3.88(l to express the integral in terms of the 
variables x an d ip instead of Xg and A r . 

10.2. Derivation 

Taking a time derivative of the expression l|3.23|l for K gives 
dK 

— = u 7 V 7 (i4T Q/3 u Q u /3 ) = \J (l K aP) u 1 u a u 15 + 2i4T Q ^u Q M 7 V 7 u /3 

= 2K a pu a a p . (10.6) 

Here we have used the Killing tensor equation V^K a ^ — 0. Now substituting the 
expression l|2.14|l for the self-acceleration gives 

(XT /I 

= A— 1 K W a V^'Prad- (10-7 

fj, dr fj, 

Now the first term here is a total time derivative to leading order in q, since K is 
conserved to zeroth order in q. Hence this term can be neglected for the reason 
explained in Sec. 19.31 Dropping the first term and substituting into Eq. (|9.6|) we get 

(f ) t = |(^VV,<W> A . (10.8) 
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Next, we use the explicit expression (19.16(1 for the radiation field ^rad- This gives 



\ dt 1 1 16nifiT ^ |w mfcn | 



E 



^mkn 



1 



| "A 
2 



^mknPmkn \t~a\ 7 down /vz^Q/3„, t-7 down\ 
\UmknPmkn\ |PAf 



We now define the amplitudes 



out 



and 



yd own 
Imkn 



2nq 



(10.9) 

(10.10) 
(10.11) 



Substituting these definitions into Eq. 110.9(1 gives the formula ((10. 1|) . Therefore it 
remains only to derive the formula 110.2(1 for the amplitudes Z? ut and Z^ own . We will 
derive the formula for Z^; the derivation of the formula for Z^ own is similar. 

We start by simplifying the differential operator K af3 u a V p which appears in the 
definitions 1(10.10(1 and ((10.11(1 . Using the expression ((3.25(1 for the Killing tensor we 
can write this as 



K als u a Va = E(r«„)/V p + S(n Q u Q )/ /3 V /3 + r 2 — . 

ut 



(10.12) 



Using the definitions 1(3.7(1 and ((3.8(1 of I and ft, the definitions ((3.14(1 and 1(3.15(1 of 

E and L z , and the fact that dt and d,p reduce to ioj m kn and — im when acting on 
(*CLz m )* gives 



K af} u a V 



a V p 



r~ -r- E + aL z + u r A)(— w <jj m kn + om — iAd r ) 



dr 2 A 



d i 

= r 2 - —(-w 2 E + aL z )(-m 2 uj mkn + am) - Au r d r . 

dr A 

Consider now the contribution of the first term in Eq. ((10.14(1 to the expression 
((10.10(1 for Z°^ kn . Using the relation 1(3.21(1 between dr and dX we can write this as 



(10.13) 
(10.14) 



2irq 



Er — (7T. 
dr 



1 \ = / r 2 — (Tr° ut I 



2?rg 

~r7 



dA 



(10.15) 



We can integrate this by parts with respect to A; the boundary term which is generated 
can be neglected for the reason explained in Sec. 19.31 This gives 



2nq I dr 



Ti \ d\ 



-fa 



uj m k n lm 



2nq 



(2iru r A(^ kn i m T) x , (10.16) 



where we have used dr/dX = £u r = Au r from Eqs. 1(3.21(1 and ((3.1|) . Combining this 
with the result obtained from substituting the second and third terms in Eq. ((10.14(1 
into Eq. ((10.10(1 gives 
2nq 

~F 

where the functions G m k n (X r , Xg) and G(X r , Xg) are [cf. Eqs. ((10.4(1 and 1(10.5(1 above] 
G m kn(X r , Xg) = - — {vj 2 E - aL z )(w 2 uj m kn - am) + 2iru r A, (10.18) 



Z 



h„k„ - — {G,nkn(X, X)(n° u ^ knlm )* + G(A, A)5 r (7r°" t fcn/m )*) A , 



(10.17) 
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G(X r , Xg) = — iAEu r . (10.19) 

As explained above, it is understood that the functions £(r, 9), A(r) and w{r) on the 
right hand sides of Eqs. (|10.18fl and (|10.19|l are evaluated at r = r(X r ) and 9 — 9(Xg). 
Also it is understood that u r is given as a function of r by Eq. 13.27( 1. and hence as 
a function of A r by using r = r(X r ) and by resolving the sign ambiguity in Eq. (|3.27() 
using Eq. 

The average over A in Eq. 110.17|l can be evaluated using the same techniques as 
in Sees. |$1 and [5] above. Using the definition l|5.43() of the "out" mode function and 
the definition (|8.9|l of uj m k n , the quantity inside the angular brackets in Eq. (|10.17fl 
can be written as 

J mk n(X,X), (10.20) 

where 

j mk n(X r ,Xe) = 2e~ m ^e^ t ° e^ kr ^ +nr ^Q Urnknlm [9(X )}* 

x e -irnA4> r (\ r ) e -imAcf) e (\ e ) e iu mkn At r (\ r ) e iuj mkn At 9 (\ d ) 

{^out,down ^ 
G mkn [X r , \e]R ZtZMK)V + G[X r , X 8 ] [r(X r )]* J . 

(10.21) 

This function is biperiodic: 

Jmkn(X r + A r , Xg) — J m kn{X r , Xg) 

Jmkn(X r , Xg + Ag) = J m kn{X r , Xg) . (10.22) 

Hence it can be expanded as a double Fourier series, and the average over A of 

Jmkn(X, A) is 

(j m kn(X,X)\ = — — [ dX r [ dXg J mkn (X r , Xg) . (10.23) 
\ / A A r Ag J Q J 

Substituting Eq. (|10.21|) into Eq. !|10.23fl and back into Eq. I|10.17|) now gives the 
formula (JTU3. 

11. Circular orbits 

11.1. Known formula for time derivative of the Carter constant 

In this section we study the prediction of the formula 1)10. l|l for the time derivative of 
the Carter constant, for the special case of circular orbits (orbits of constant Boyer- 
Lindquist radius r). It is known that circular orbits remain circular while evolving 
under the influence of radiation reaction in the adiabatic regime |52 EH E3|. Also 
for circular orbits it is possible to relate the Carter constant K to the energy E and 
angular momentum L z of the orbit: specializing Eq. I|3.26fl to zero radial velocity 
u r = yields 

K = ~ 1 -^--r 2 . (11.1) 

In order for this relation to be preserved under slow evolution of K, E and L Zl we 
must have 

'dK\ _2(m 2 E-aL z ) \ 2 /dE\ / dL z s 



dt l t 



dt 1 1 \ dt 



(11.2) 



We now verify that our expression (|10.1|) for (dK/dt) t satisfies the condition 111.211 . 
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11.2. Verification that our result reproduces this formula 

For circular orbits, the following simplifications apply to the harmonic decomposition 
of the geodesic source derived in Sec. |S] First, the functions Ai r (A), t r (X), A</> r (A) 
and r (A) defined by Eqs. (G3EJ), EHI, COSt and iPOk^I vanish identically. Second, 
the biperiodic function J u im(X r ,Xe) defined in Eq. i|8.17f) is therefore independent of 
A r . It follows that the double Fourier series i|8.2U[) is replaced by the single Fourier 
series 

oo 

Jui m (Xe)= J^rn k e- lkTl>Xl> , (11.3) 

k— — oo 

where the coefficients J^imk are given by 

Julmk = -J- / " d\ e e lkT ° x °,U m (Xe). (11.4) 
A e Jo 

The effect of this change on the subsequent formulae for fluxes and amplitudes in Sees. 
l^land lTUl can be summarized as: (i) remove the sums over n; (ii) remove the indices n 
from the amplitudes Zi m j. n ; (iii) remove the averaging operation 

1 f K 
— \ dX r 

A r Jo 

from the formulae for the amplitudes; and (iv) evaluate the integrands in the formulae 
for the amplitudes at n — 0. Also the formula l|8.9|) for the frequency oj m kn is replaced 

by 

W m fc = mSlcj, + kVLg. (H-5) 

With these simplifications, the formula (|9.7(l for the time derivatives of energy 
and angular momentum becomes 



^ mk 



d£ \ _ 1 \ -« f U m k 

dt/t ~ 64tt 2 ^ \ m J \ Umk \ 



1 | ryOMt | 2 



I |2 I Irak 

\0>lmk\ 



^mkPmk \ T lmk\ i^down|2 



2 I link 



(11.6) 



where £ — E or L z , and the expression in curly brackets means either cj m fc (for _E) 
or to (for L z ). Similarly the expression l|lU.lfl for the time derivative of the Carter 
constant becomes 

' / yout \ * yout 

■K^lmk) ^Imk 



\dt/ t ~ 327r 2 //^Kfc| 



|a; m fc| 2 

2 



^mkPrnk IVmkl / ydown \ * ydown /-i -i -rN 

+ I \Ta HV^lmk ) ^Imk • \ Ll -<) 

\<jJmkPmk\ \Plmk\ 

Comparing Eqs. I|11.6[) and l|11.7fl . we see that the condition l|11.2|l will be satisfied as 
long as the amplitudes Zi m k satisfy 



70ut,down V 137 E a Lz) 1 2 \ rzout.down 

hmk = ^ \ W U ™k - ma) Z lmk 



Therefore it suffices to derive the identity l|11.8[) . We will derive this identity for the 
"out" modes; the "down" derivation is identical. 



Computing inspirals in Kerr in the adiabatic regime 



42 



The particular form of the expressions for the amplitudes that are most useful 
here are Eqs. (|9~22|) and 110.170 : 

« = -y(s(Cl)*) A , (11.9) 

and 

ZZ\ = ^ (G mk (X, A) (ttZIY + G(\, X)d r K U A-)*) A . (11.10) 

Here we have defined = 7T° ut lm . It follows from Eq. I|10.19|) together with u r = 
that the function G(A, A) vanishes, so the formula for Z°^ k simplifies to 

ZZl = ^ (Gmk(\ A) (nZ t kY) x • (11-11) 

We now substitute into Eq. (|ll.llf> the expression (|10.18fl for G m k- The second term 
in Eq. (|10.18fl vanishes since u r — 0, and of all the factors in the first term, only £ 
depends on A; the remaining factors are constants and can be pulled out of the average 
over A. This yields 

- t _ 2ng {w 2 E-aL z ) „ m WvLout\A m 19^ 

Z lmk f ^ \ mk ~ I \ y^lmk) / • (.H- 12 ) 

Comparing this with the formula l|11.9fl for Z°£\ now yields the identity l|11.8fl . 
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Appendix A. Accuracy of adiabatic waveforms 

In this appendix we estimate the accuracy of the adiabatic waveforms, expanding 
on the brief treatment given in Ref. 32 . Specifically we estimate the magnitude of 
the correction to the waveform's phase that corresponds to the term 4>2{t) in Eq. 
I|1.2|l . We can roughly estimate this term by using post-Newtonian expressions for 
the waveform in which terms corresponding to ^(t) are readily identified. While the 
post-Newtonian approximation is not strictly valid in the highly relativistic regime 
near the horizon of interest here, it suffices to give some indication of the accuracy 
of the adiabatic waveforms. In this appendix we specialize for simplicity to circular, 
equatorial orbits. We also specialize to gravitational radiation reaction, unlike the 
body of the paper which dealt with the scalar case. Our conclusion is that adiabatic 
waveforms will likely be sufficiently accurate for signal detection. 

The procedure we use is as follows. We focus attention on the phase ^(f) of the 
Fourier transform of the dominant, I = 2 piece of the gravitational waveform. Here 
/ is the gravitational wave frequency. Now a change to the phase function of the 
form ^(Z) — > \&(/) + + 27r/At is not observable; the constant At corresponds to 
a change in the time of arrival of the signal. Therefore we focus on the observable 
quantity d 2 ^/df 2 , and we write this as 

^l(f) = G(f, mi ,m 2 ). (A.l) 
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where 



and 



Here G is a function, m\ is the mass of the black hole, and m 2 is the mass of the 
inspiralling compact object*. The specific form of G obtained from post-Newtonian 
theory is discussed below. 

We next expand the function G as a power series in m 2 , at fixed / and mi. The 
leading order term scales as m 2 l , and we obtain 

G(f, m u m 2 ) = G ° (/ ' mi) + G 1 (f, m x ) + G 2 (f, mi )m 2 + . . . (A.2) 
m 2 

The first term in this expression corresponds to the term &i(t) in Eq. Q1.2|) : it gives the 
leading-order, adiabatic waveform. The error incurred from using adiabatic waveforms 
is therefore 

AG(f, mi ,m 2 ) = G(J, mi ,m 2 ) - G °^' m ^ = Gl (/ ; TOl ) + G ^f mi ) TO2 + . . . ( A .3) 

m 2 

We want to estimate the effects of this phase error. 
It is useful to split the phase error into two terms: 

AG(/,m 1; m 2 ) = AG 1 (f,m 1 ,m 2 ) + AG 2 (f,m u m 2 ), (A.4) 

An if \ nit \ Go(/,TOi+m 2 ) 

AGi(/,mi,m 2 ) = G(/,mi,ra 2 ) (A. 5) 

m 2 

AG 2 (f,m 1 ,m 2 ) = . (A. 6) 

m 2 m 2 

The phase error AG 2 corresponds to the error that would be caused by using an 
incorrect value of the mass of the black hole. The effect of this error on the data 
analysis would be to cause a systematic error in the inferred best fit value of the black 
hole mass. However, the fractional error would be of order fi/M <§; 1. This error will 
not have any effect on the ability of adiabatic waveforms to detect signals when used 
as search templates. Therefore, we will neglect this error, and focus on the remaining 
error term AGi. 

We now compute explicitly the phase error A^i(f) that corresponds to the term 
AGi in Eq. i|A.4|) . We use the post-3. 5-Newtonian expression for ^(J), which can be 
computed from the orbital energy E(f) and gravitational wave luminosity E{f) via 
the equation j^S] 

f = -™ 

For the energy E(f) we use the expression given in Eq. (50) of Ref. [SSjj and for the 
luminosity E we use Eq. (12.9) of Ref. |H7] j). This gives 

3M 
128/V 

* We also define M = m\ + m,2, the total mass, and fi = mim2/M, the reduced mass. Throughout 
the body of the paper we referred to the mass of the inspiralling object as fi and the mass of the 
black hole as M; this is valid to leading order in fi/M. In this appendix however we need to be more 
accurate, and hence we revert to using mi and m,2 for the two masses. 

(t With parameter values A = -1987/3080 from Ref. gHl and 6 = -11831/9240 from Refs. |59ll60| . 



*(/) = 2^/< c + 4>c + 777ET75^(/)> ( A " 
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where t c and 4> c are constants, y = (ttM/) 1 / 3 , and 

3715 y 2 o 15293365 y 4 11583231236531 y 6 6848 7 y 6 

■F = H 16 Try 1 



756 " 508032 4694215680 21 

640 7r 2 y 6 77096675 Try 7 55y 2 /z 27145 yV 

3 + 254016 + 9M + 504 M 
15737765635 y 6 n 2255 tt 2 y 6 /i 1014115 Try V 3085 y 4 ^ 2 

3048192 A/ + 12 M + 3024 M + 72 M 2 
76055 yV 2 36865 iry 7 n 2 127825 yV ay 3 / 2 75/z 
+ 1728 M 2 378 Af 2 1296 A/ 3 +_ 3\ 11 ^ + ~aT 



13696 y 6 log(2) , 38645 tt y 5 log(y) 6848 y 6 log(y) 

(A.9) 



21 252 21 

5 7ry 5 /i log(y) 



M 

Here 7 is Euler's constant, a is the dimensionless spin parameter of the black hole, and 
v = mi/M = 1/2 + ^/1/4 - fi/M. We have also added to Eq. (|A.9|I the spin-dependent 
terms up to post-2-Newtonian order, taken from Eq. (1) of Ref. [HI], specialized to a 
non-spinning particle on a circular, equatorial orbit. 

We now insert the expression (|A.8|I for ^(f) into Eqs. i|A.l|) and (|A.5|I . and 
integrate twice with respect to / to obtain A^i(f). The result is 

"55 y 2 27145 y 4 15737765635 y 6 2255 7r 2 y 6 1014115 Try 7 



128y 5 



9 504 3048192 12 3024 

3085 yV 76055 y 6 /i 36865 Try 7 ^ 127825 yV 

72 M + 1728 M 378 M 1296 M 2 

3 



-5 Try 5 log(y) + ^-(75 + 113 g ) 



(A.10) 



3 

where q — (y 2 — l)M/(i = — 2 + 0(fj,/M). Now as discussed above, a change to the 
phase function of the form ^(f) — > + ^0 + 2^/ At is not observable. This freedom 
allows us to set *(/o) = ^'(/o) = where /o is a fixed frequency, which amounts to 
replacing our expression A^i(/) for the phase error with 

S^(f) = A*i(/) - A^xCfo) - A*i(/o)(/ - /o). (A.11) 

We now evaluate the phase error for some typical sources. First, we take 

/j, = 1M Q , M = lO 6 A/ , a = 0.999. For this case the last year of inspiral extends 
from / = 10~ 2 Hz to ~ 3 x 10 -2 Hz ^21; the maximum value of \S^i\ is 0.19 cycles 
if we take /q = 0.018 Hz. As a second example with a higher mass ratio, we take 
/Lt = IOMq, M = 10 6 Af Q , a = 0.999. The last year of inspiral for this case extends 
from / = 0.004 Hz to / = 0.03 Hz [T2]; the maximum value of |<5*i| is 0.47 cycles 
if wc take f = 0.0138 Hz (see Fig.ETJtt- For signals fr om intermediate mass black 
holes that may be detectable by LIGO, the phase errors are yet smaller. 

These phase errors are sufficiently large that adiabatic waveforms cannot be used 
as data-analysis templates; the template will not quite stay in phase with the signal 
over the entire ~ M/jj, ~ 10 5 cycles of inspiral. However, the requirements on detection 
templates are much less stringent: phase coherence is only needed for ~ 3 weeks, rather 
than a year 0]. In addition, in the matched filtering search process phase error will 
tend to be compensated for by small systematic errors in the best-fit mass parameters. 

ff If the spin-dependent term in Eq. IA.10I is omitted, the errors are about ~ 50% larger. 
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Phase error from adiabatic waveforms 

0.2 r-, , , , , 1 



0.1 - 




0.005 0.01 0.015 0.02 0.025 0.03 
f (gravitational wave frequency, Hz) 

Figure Al. This plot shows the error in the phase of the Fourier transform of 
the waveform that results from using adiabatic waveforms, estimated using post- 
3.5-Newtonian waveforms, for a lOM© compact object spiralling into a 10 6 Mq 
black hole with spin parameter a = 0.999. For this system the last year of inspiral 
extends from / = 0.004 Hz to / = 0.03 Hz |T2|. 



Because the adiabatic waveforms are almost good enough for data-analysis templates 
(phase errors ~ 1 cycle in the worst cases) , adiabatic waveforms will likely be accurate 
enough for detection templates both for ground-based and space-based detectors. 

Appendix B. Typos in Gal'tsov 

This appendix lists some of the typos in the paper by Gal'tsov |36) : 

• The discussion of the range of summation of the indices I, m given after Eq. (2.5) 
is incorrect. It should be summation l,m with I > \s\ and \m\ < I rather than 
/ < \s\ and |m| < 1. 

• In the definition of the variable s u given just before Eq. (2.7), the factor 
(r 2 + a 2 )- 1 ' 2 should be (r 2 + a 2 )' 1 / 2 . 

• In the first of the four equations in Eqs. (2.9), the term r 2 e~ tur should be 

• The right hand side of Eq. (2.25) should be divided by — 47r, as already discussed 
in Sec. 16. II above. 

• In Eq. (2.30), the quantity f? ut should be replaced by its complex conjugate f? ut . 

• In the first line of Eq. (3.2), the quantity Z(6',ip') should be replaced by its 
complex conjugate Z($',<p')*. 

• In the second of Eqs. (3.3), the right hand side should read 

K s T s ? s (wk/\wk\)( s v up + a- s s v in ) 

rather than 

K s T s ? s (wk/\wk\)( s v up ) + a- s s v in , 
i.e. the closing bracket is in the wrong place. 
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• In the second term on the first line of Eq. (3.6), the first factor of s v up (r) should 
be omitted. 

• In Eq. (3.9), the argument of the last factor should be x' rather than xx' . 

• Gal'tsov's equation (4.13) for the rate of change of energy or angular momentum 
of the orbit is correct only for sources which are smooth functions of frequency. 
However, for bound geodesies, the inner products (n, J) are sums of delta 
functions in frequency, so it does not make sense to square these inner products 
as Gal'tsov does. The correct version of this equation is our Eq. I|9.7|l above. 
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